University of Science and Technology of China 

Hefei, CHINA 



Measurement of charm production cross-section and 
leptons from its semileptonic decay at RHIC 



A dissertation submitted 
for the degree of 
Doctor of Philosophy in Physics 

by 

Yifei Zhang 



Co-supervisors: Ziping Zhang, Nu Xu 



2007 



(c) Copyright by 
Yifei Zhang 
2007 

All Rights Reserved 



1 



Dedicated to my dear family 



2 



Acknowledgments 



I would never forget people whom I was/am working with. And I can not finish my 
Ph.D project without their support and help. This thesis includes many efforts and 
contributions from them. 

Firstly, I would thank Prof. Ziping Zhang and Prof. Hongfang Chen for guiding 
me into the heavy ion collision physics field and offering me the great opportunity to 
participate the STAR data analysis. I appreciate their continuous supervision and support 
in the last six years. I thank Dr. Hans-Georg Ritter and Dr. Nu Xu who provide me the 
environment and opportunity to work with the great RNC group in LBNL. I thank Dr. 
Nu Xu for his guidance and fruitful discussions on many aspects of my research in the last 
two years. I greatly benefited from Dr. Nu Xu's profound and broad understanding of 
physics. I am deeply inspired by his assuredness and encouragement. I would especially 
thank Dr. Zhangbu Xu from BNL for his rich experience in data analysis and patience in 
instructing me all the analysis details throughout the thesis. 

I would like to thank Dr. Haibin Zhang and Dr. Chen Zhong for their helpful coop- 
eration in BNL and plentiful discussions. Thanks also go to Dr. Xin Dong, Dr. Haidong 
Liu, Dr. Yan Lu and Dr. Lijuan Ruan for their amount of discussions on work and for 
the joyful time we had together. 

I am very lucky to work with the RNC group in LBNL. I am impressed by the folksy 
and enthusiastic people in the group and their open mind and meticulosity on research. I 
learned a lot from the fruitful discussion with them. Especial thanks to Dr. Jim Thomas, 
Dr. Ernst Sichtermann and Dr. Andrew Rose who gave me lots of help and useful 
discussions on the HFT software. I also appreciate Dr. Marco van Leeuwen and Dr. Vi 
Nham Tram for their valuable discussion on charm physics and help on the correction of 
charm paper and my thesis. I also thank other members in RNC group for their kind 
assistances. 

I would thank all the members of spectra and heavy flavor working group. Especial 
thanks to Dr. Bedanga Mohanty for the plenty of discussions on many aspects of physics. 



3 



I would thank to Prof. Huan.Z. Huang, Dr. Xiaoyan Lin from UCLA for many helpful 
discussions on heavy flavor physics. I also thank Dr. Thomas Ullrich, Dr. Manuel, Dr. 
Weijiang Dong and other people who gave me lots of useful suggestions and discussions 
on the detailed data analysis of heavy flavor. 

I would especially thank Prof. Charles Whitten from UCLA for his elegant wording 
and grammatical corrections on my thesis. 

Many thanks to our high energy physics group (HEPG) in USTC. Especial thanks 
to Prof. Xiaohan Wang, Prof. Zizong Xu, Prof. Jian Wu and Dr. Ming Shao for their 
continuous discussions and support. I would thank Dr. Yane Zhao and Dr. Yi Zhou for 
their hospitable help on my living. And many thanks to my classmate Dr. Qing Shan, 
Zebo Tang and other students in our group for their hardwork on the MRPC production 
and testing. 

Thank the PDSF and RCF software group for their offered resources and continuous 
development and support on software. Thank the whole STAR Collaboration for building 
such a beautiful detector and provide the enormous amount of data with high quality. I 
would particularly thank STAR TOP group and HFT group for the hardwork on detector 
upgrades and brightening the future. 

Finally, what I will remember to the end of my life is the eternal love, the toler- 
ant consideration and selfless dedication from my family. Their continuous support and 
understanding have always been my motivity of never giving up. 



4 



Abstract of the Dissertation 

Measurement of charm production cross-section and 
leptons from its semileptonic decay at RHIC 

by 

Yifei Zhang 

Doctor of Philosophy in Physics 
University of Science and Technology of China, Hefei, 2007 

The strong interaction, one of the four fundamental forces of nature, is confining: there 
is no single quark as a color-triplet state observed experimentally. Quantum Chromody- 
namics (QCD) is believed to be a correct basic gauge field theory of strong interactions. 
Lattice QCD calculations predict a phase transition from hadronic gas to a new matter. 
Quark- Gluon Plasma (QGP), at high energy nuclear-nuclear collisions. In this new form 
of matter, quarks are deconfined and approach local thermalization. One of the ultimate 
goals of the heavy ion collision experiments is to search for the QGP matter and study its 
properties. The Relativistic Heavy Ion Collider (RHIC) located at Brookhaven National 
Laboratory (BNL) provides a high energy density environment to create and search for 
the QGP matter by colliding ions hke Au at energies up to y^s^=200 GeV. Recent ex- 
perimental studies at RHIC have given strong evidences that the nuclear matter created 
in Au-|-Au collisions at y/s^=200 GeV has surprisingly large collectivity and opacity as 
reflected by its hydrodynamic behavior at low px and its particle suppression behavior at 
high pt. 

Charm quarks provide a unique tool to probe the partonic matter created in relativistic 
heavy- ion coUisions at RHIC energies. Due to their large quark mass (~ 1.3 GeV/c^), 
charm quarks are predicted to lose less energy than light quarks via only gluon radiation. 
A measurement of the nuclear modification factor for the charmed hadrons semileptonic 
decayed single electrons compared to light hadrons is valuably important to complete 
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the picture of the observed jet-quenching phenomenon and help us better understand the 
energy-loss mechanisms at parton stage in Au-|-Au collisions at RHIC. 

Furthermore, the interactions between charm quarks and the medium could boost 
the radial and elliptic flow resulting in a different charm px spectrum shape. Due to its 
large mass, a charm quark could acquire flow from the sufficient interactions with the 
surrounding partons in the dense medium. The measurement of charm flow and freeze- 
out properties is vital to test light flavor thermalization and the partonic density in the 
early stage of heavy ion collisions. 

Charm quarks are believed to be produced only at early stages via initial gluon fusions 
and its cross section can be evaluated by perturbative QCD calculations. Thus study of the 
binary collision (A'bin) scaling properties for the total charm cross-section from different 
collision systems can test the theoretical assumptions and determine if charm quarks are 
indeed good probes to the partonic matter created in high energy heavy ion collisions. 
Charm total cross-section measurement is also essential for the separation of bottom 
contribution in non-photonic electron measurement, and for the model calculations, which 
tries to explain the observed similar suppression pattern of J/^ at RHIC and SPS. 

In this thesis, we present the measurements of Kn at low pr {pt < 2 GeV/c) 

and non-photonic electron spectra (0.9 < pr ^ 5 GeV/c) from semi-leptonic decay. 
In addition, we use a newly proposed technique to identify muons from charm decays at 
low Pt- The combination of all these three measurements stringently constrains the total 
charm production cross-section at mid-rapidity at RHIC. They also allow the extraction 
of the charmed hadron spectral shape and a study of possible charm radial flow in Au+Au 
collisions. 

D° mesons were reconstructed from hadronic decay D° — > K~7r'^ — > X+tt") with 
a branching ratio of 3.83% in minbias Au-|-Au collisions. The yields were obtained 
from fltting a Gaussian plus a linear (or a second-order polynomial function) for resid- 
ual background to the invariant mass distributions of kaon-pion pairs after mixed-event 
combinatorial background subtraction. 
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The inclusive muons at 0.17 < Pt < 0.25 GeV/c in minbias Au+Au and central 
Au+Au collisions were analyzed by combining the ionization energy loss (dE/dx) mea- 
sured in the Time Projection Chamber (TPC) and the mass calculated from the Time Of 
Flight (TOF) detector at STAR after the residual pion contamination subtracted. The 
dominant background muons from pion/kaon weak decays were statistically subtracted 
using the distribution of the distance of closest approach to the collision vertex (DCA). 

Inclusive electrons up to pT = 5 GeV/ c are identified by using a combination of velocity 
measurements from the TOF detector and dE/dx measured in the TPC. Photonic 
background electrons are subtracted statistically by reconstructing the invariant mass of 
the tagged and every other partner candidate using a 2-dimensional invariant mass 
method. Partner track finding efficiency is estimated from STAR Geant + Monte Carlo 
embedding data. More than ~ 95% of photonic background (photon conversion and tt" 
Dalitz decay) can be subtracted through this method. The pt spectra for non-photonic 
electrons in central Au-|-Au minbias Au-|-Au collisions and its subdivided centralities will 
be presented. In addition, the method to extract the elliptic flow V2 of non-photonic 
electron is also developed in the thesis. 

By combining these three independent measurements: D — > Kn , muons and elec- 
trons from charm semileptonic decays in minbias and central Au-|-Au collisions at RHIC, 
we observed that: The transverse momentum spectra from non-photonic electrons are 
strongly suppressed at 0.9 < Pt < 5 GeV/c in central Au-|-Au collisions relative to d+Au 
collisions. For electrons with pr ^ GeV/c, corresponding to charmed hadrons with 
Pt ^ 4 GeV/c, the suppression is similar to that of light baryons and mesons. The blast 
wave fit to the electron spectra with Pt ^ 2 GeV/c indicates that charmed hadrons may 
interact and decouple from the system differently from multi-strange hadrons and light 
hadrons. Future upgrades with a direct reconstruction of charmed hadrons are crucial for 
more quantitative answers. Charm differential cross-sections at mid-rapidity {da^^ /dy) 
are extracted from a combination of the three measurements covering 90% of the kine- 
matics. The total charm cross-section per nucleon-nucleon collision (c^^) is reported as 
lA0±0.11{stat.) ±0.39{sys.) mb in 0-12% central Au+Au and 1.29 ±0.12 ±0.36 mb in 
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minbias Au+Au collisions at y/s^=200 GeV. The charm production cross sections are 
found to follow the number of binary collisions scaling. This supports the assumption that 
hard processes scale with binary interactions among initial nucleons and charm quarks 
can be used as a probe sensitive to the early dynamical stage of the system. 

In the above measurements, the bottom contribution to the non-photonic electron 
spectrum is neglected. The separation of bottom and charm contributions in current 
non-photonic electron measurements is very difficult. There are large uncertainties in 
the model predictions for charm and bottom production in high-energy nuclear collisions. 
Thus identification of bottom from the non-photonic electron measurements is crucial to 
better understand charm physics. In the discussion section of this thesis, we will try a fit 
to non-photonic electron spectrum and estimate the bottom contributions. We will also 
compare the V2 distribution from simulation to the experimental data and estimate the 
possible charm V2- 
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CHAPTER 1 



Introductions 

1.1 Quantum ChromoDynamics 

1.1.1 Confinement and asymptotic freedom 

Quantum ChromoDynamics (QCD) |DKS03j is tiiought to be a correct theory of the 
strong nuclear force, one of the four fundamental forces of nature. It describe the strong 
interactions among quarks, which are regarded as fundamental constituents of matter, 
via their color quantum numbers. The strong interactions among quarks are mediated 
by a set of force particles known as gluons. Different from Quantum ElectroDynamics 
(QED) - the gauge theory describing electromagnetic interaction, QCD is based on the 
non-Abelian gauge group SU{3), with gauge bosons (color octet gluons), and hence the 
gluons could have self-interacting. This results in a negative /^-function and asymptotic 
freedom at high energies and strong interactions at low energies. 

These strong interactions are confining: the self-coupled gluons strongly restrain the 
isolation of the quarks at large distance. There is no single quark as a color-triplet state 
observed experimentally. Only color-singlet bound states can propagate over macroscopic 
distances. The only stable color-singlets with size of the order of Ifm are quark— antiquark 
pairs, mesons, and three-quark states, baryons. At high energy reactions, like deep in- 
elastic scattering, the quark and gluon constituents of hadrons act as quasi-free particles, 
partons. Such reactions can be factorized into the convolution of non-perturbative parton 
distribution functions, which cannot be calculated from first principles directly. But with 
process-dependent functions {i.e. hard processes involving large momentum transfers), 
the reactions can be calculated as perturbative expansions in the coupling constant as- 
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In QED, the electrodynamic coupling constant a = j^. However, due to the gluons 
self-interactions, the renormalized QCD coupling shows renormalization scale (/i) depen- 
dence |Alp02| . The running coupling q;s(/x) can be written as: 



where gs, which is strong charge in the gauge group, is the only parameter in the QCD 
Lagrangian besides the quark masses. (3q (>0) is the first coefficient of the 6eta-function 
(renormalization neglects the higher orders). The strong force of the gluon-gluon self- 
coupling becomes smaller at shorter distance or with larger momentum transfers (a^^ 
as yU — > oo), which is known as asymptotic freedom. In this case, QCD can be calculated 
perturbatively. Many experiments measured at different scales. Since some of the 
precise measurements come from decays, it has become universal to use as{Mz) as 
the label. The as{Mz) = 0.1176 ± 0.002 |PDG04j comes from a fit to the experimental 
data, and the QCD scale Aqcd ~ 200 MeV. Fig. 11.11 shows the measured as at different 
momentum transfer scale fi compared with Lattice QCD calculations. 

1.1.2 Nonperturbative QCD 

Corresponding to the asymptotic freedom of QCD at high momentum scales (short dis- 
tances), the running coupling becomes larger at lower momenta (longer distances). This 
effect is called infrared slavery. Eg. 11.11 shows strong coupling at /i ~ Aqcd- Perturbation 
theory is not expected to be a reliable solution when as becomes large. Thus nonper- 
turbative methods must be used to study the low-momentum, long-distance behavior of 



A closely related phenomenon is hadronization. This term refers to the dynamical 
process in which a set of partons produced in a short-distance interaction reorganize 
themselves, possibly with the production of additional partons, to make the hadrons that 
are observed in the final state. This again is expected to be a low momentum transfer 
process, corresponding to the relatively long timescale over which hadron formation takes 
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Figure 1.1: Measured QCD running coupling constant from different experiments 

compared with Lattice QCD calculations. 
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A nonperturbative approach that is widely used to study confinement is lattice QCD |Gup98| . 
The continuum field theory is replaced by one defined on a space-time lattice. As long as 
the size of the lattice elements inside a hadron is sufficiently small, hadronic properties 
and matrix elements should be reliably calculable on massively parallel computers. 



1.1.3 Heavy quark masses and perturbative QCD (pQCD) 

The masses of the heavy c quark and b quark {iric ~ 1.3 GeV/c^, mf, ~ 4.75 GeV/c^) 
are almost exclusively generated through their coupling to the Higgs field in the electro- 
weak sector, while the masses of light quarks (u, d, s) are dominated by spontaneous 
breaking of chiral symmetry in QCD. This means that in a Quark-Gluon Plasma (QGP), 
where chiral symmetry is expected to be restored, light quarks are left with their bare 
current masses while heavy-flavor quarks remain heavy. The heavy quark masses are not 
modified by the QCD vacuum jPDG04[ISe04j . Hence once heavy quarks are pair-produced 
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in early stage via initial gluon fusions [Lin95j at high energy collision, their total number 
is conserved because the typical temperature of the subsequently evolutional medium is 
much smaller than the heavy quark masses. Therefore, heavy-flavor quarks are an ideal 
probe to study early dynamics in high-energy nuclear collisions, and their cross section 
can be evaluated by pQCD calculations. 




u d s c t b 

Figure 1.2: Masses of the six flavors. The masses generated by electroweak symmetry 
breaking (current quark masses) are shown as dark boxes; A large fraction of the light 
quark masses is due to the chiral symmetry breaking in the QCD vacuum (constituent 
quark masses), shown as light boxes. The numerical values were taken from 

Ref. |PDG04j . 

QCD is formulated in terms of quarks and gluons while the experimentally observed 
states are hadrons. As previously mentioned, in principal the method of perturbation 
theory is appropriate in the high-momentum scale, short-distance regime. It is known 
as the QCD factorization theorem |Css89t ICol98] and states that for processes that have 
initial and/or observed final state hadrons the differential cross-section has the following 
form: 

da{x, Q^, m^) = Y[Y^ fi/h{x, fi^) ® dai^f{x, Q^ m^ fi^, /i|) ® Dh,/f{x, fi^) + 0{A/Q), 

h,h' i,f 

(1.2) 

where the factor fi/^ stands for the Parton Distribution Function (PDF) of the parton i 
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inside the hadron h present in the initial state, and Q and x represent the hard scale and 
some kinematical variable respectively. The parton distributions depend on the factor- 
ization/renormalization scale fi"^. The second factor dai^f, also known as the (Wilson) 
coefficient function, represents the partonic hard scattering cross section for the reaction 
i ^ f that depends on the unphysical renormalization and factorization scales /i^ and fi'jp 
and on the masses of the heavy quarks m^. The last factor D is the so called Fragmenta- 
tion Function (FF). It contains the information for the hadronization of the hard parton / 
(that is produced in the hard process described by the partonic cross-section da) into an 
observed hadron h'. Since the concepts of parton distributions and hadronization are only 
for the initial and final states of hadrons, lepton-nucleon Deep Inelastic Scattering (DIS) 
and high energy e^e^ collisions are performed to measure PDFs and FFs, respectively. 

The large masses (^ Aqcd) of the heavy quarks, which are 'hard' scales, make per- 
turbative QCD applicable | Vog02| . Asymptotic freedom is assumed in calculating the 
interactions between two heavy hadrons on the quark/gluon level but the confinement 
scale determines the probability of finding the interacting parton in the initial hadron. 
Factorization assumes that between the perturbative hard part and the full expansion, 
non-perturbative parton distribution functions are applied. The integral hadronic cross 
section in an AB collision where AB = pp, pA or nucleus-nucleus can be written as: 

aAB{S,m'^)= ^ f — j dxidxj5{xiXj-T)ff'{xi,ii],)ff{xj,iJ.],)^ 

(1.3) 

where ff" and are the nonperturbative PDF determined from experiments, Xi and Xj 
are the fractional momentum of hadrons A and B carried by partons i and j , t = s/ S , s 
is partonic center of mass energy squared, (Jjj(s, m^, /x^) is hard partonic cross section, 
which only depends on quark mass m, not kinematic quantities, can be calculated in QCD 
in powers of known as leading order (LO), n = 0; next-to-leading order (NLO), n = 1 
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.... Eq. 11.41 express the partonic cross-section to NLO: 



(1.4) 



where p = Am? / s, ppis factorization scale which separates hard part from nonperturbative 
part, fiR is renormahzation scale, at which the strong coupling constant as is evaluated. 
Pf = /^_R is assumed in evaluations of parton densities, and f-j'''^ are dimensionless, 
/i— independent scaling functions. For LO, a = b = and ij = qq, gg. For NLO, a = 1 
,b = 0,1 and ij = qq,gg and qg,qg. fij'^^ are always positive, f^j'''^ can be negative. If 
p?p = m?, f-j'^^ does not contribute. Results of theoretical calculations strongly depend 
on quark mass, m, factorization scale, pp, in the parton densities and renormahzation 
scale, PR, in a„ etc. |MNR93l rVog02i IRP03j . 



1.1.4 Deconfinement, equation of state and QCD phase 

Currently, the isolated quark has never been observed. Quarks and gluons are confined 
in QCD matter as hadron gas. If the temperature and energy density of the matter are 
sufficiently high, the strong force among quarks and gluons may be greatly reduced. 
Thus quarks and gluons could be deconfined from hadrons and become disengaged with 
large distance among them in that hot and dense condition. This phenomenon is called 
deconfinement. The deconfinement indicates a phase transition from a hadron gas to a 
new deconfined form of matter - the QGP, which is formed by these dissociative quarks 
and gluons with new (color) degree of freedom (d.o.f). In the QGP phase, the broken 
chiral symmetry will be restored. The evolution of energy density and pressure with 
temperature from Lattice QCD calculations illustrate the phase transition |KLP00j . see 
Fig.O 

Large increase of degrees of freedom at Tc is seen in the rapid change in energy density 
and pressure. The critical temperature is predicted to be Tc ~ 150 — 180 MeV, and 
the energy density at the critical point is edTc) ~ 1 — 3 GeV/fm^ (~ 0.17 GeV/fm^ 
for nuclear matter) |Kar02b] . The arrows indicate the Stefan- Boltzmann (SB) limits. 
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T (MeV) T (MeV) 

Figure 1.3: The evolution of e/T'^ (a) and p/T^ (b) witli the increase of temperature T 
for 3 flavor and 2 flavor configurations. The arrows indicate the SB limit for each case. 

applicable for the systems with massless, non-interacting quarks and gluons. The similar 
deviation of the curves from the SB limits indicates that besides the effect of quark 
masses, the quasiparticles, which form the new matter, must interact with one another. 
The various heavy-ion experiments devoted to the creation and detection of new forms of 
highly excited matter also indicate that, at the reachable energy scales, the produced phase 
exhibits a strong collective behavior (see next section). All these theoretical predictions 
and experimental observations are incompatible with a weakly interacting QGP. 

Based on local 4-momenta conservation and other conserved currents (e.g. baryon 
number), 

d.T^'' = 0, d,f = 0, (1.5) 
from the Hydrodynamic model, the equations of relativistic ideal hydrodynamics are 

j^'(x) = n{x)u''{x), (1.6) 

T>"' = [eix) + p(x)]m^m" - g'^'pix), (1.7) 
where e{x) is the energy density, p{x) is the pressure, n{x) is the conserved number 
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density, and u^{x) = 'y{l,Vx,Vy,Vz) with 7 = 1/ y/l—vl—Vy—v'j is the local four velocity 
of the system. 

Eq. 11.51 contains 5 equations for 6 unknown fields e,n,p,Vx,Vy,Vz. To fully describe 
the system one needs an equation of state (EOS), which relates pressure, energy and 
baryon density. The thermal condition of the system can be fixed by pressure, p and 
temperature, T{e). The EOS, which describes the system response to changes of the 
thermal condition, is separately constructed for a relativistic ideal QGP phase (dashed 
line in Fig. 11.41) . dp/de = |, and a hadron resonances gas (dotted line), dp/de ~ |. 
By combining these two EOSs, the EOS indicates the phase transition (solid line) can 
be obtained from Maxwell construction with a bag constant B. Here B is tuned to the 
desired phase transition temperature: 5^/^ = 230 MeV gives Tdn = 0) = 164 MeV at 
vanishing net baryon density. 



net baryon n=0 fm' 



Phase transition 



Hadon resonances gas 



12 3 4 

Energy density e (GeV/fm^) 

Figure 1.4: Equations of state for ideal QGP phase, hadron resonances gas and phase 

transition at vanishing net baryon density. 
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1.2 Relativist ic Heavy Ion Collisions 

With very high energy density and temperature in a bulk system, the new QGP phase 
of matter is predicted to be formed by the deconfined quarks and gluons. In order to search 
for the production of QGP matter and to study its properties, relativistic heavy-ion col- 
hsions have been performed to provide the hot and dense environment under laboratory 
conditions. The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Labo- 
ratory is built to create and search for this novel form of matter by colliding Au ions at 
energies up to >/s^=200 GeV. 




In relativistic heavy-ion collisions, due to the Lorentz contraction in the moving direc- 
tion, the collision geometry between the two nuclei behaves as two thin disks approaching 
each other nearly with speed of light. The heavy- ion coUisions can be approximated as 
the interpenetrating collisions of their constituent nucleons with partonic interactions at 
high energy. The number of participating nucleons and the produced particle multiplici- 
ties in the final hadron state can be calculated by the Glauber model which relates these 
quantities to the size of the impact parameter, b, which is defined as the distance between 
the perpendicular bisectors along the colliding direction of the two ions. The impact 
parameter is large for peripheral collisions, consisting of a small number of participants 
and small multiplicities. A central collision occurs when the impact parameter is small. 
Only central collisions of two heavy nuclei at high energies are expected to provide a suf- 
ficiently hot and dense environment to produce the QGP. In central heavy-ion collisions, 
with sufficient interactions of the partons, the chemical and local thermal equilibrium of 
the quark-gluon system could be reached and thus the QGP forms. At the initial stage 
of the collision with the primary hard scattering, high pt jets, heavy quark pairs, direct 
photons, etc. are created due to large momentum transfers. Due to high energy densi- 
ties and pressure gradients, the QCD system expands and cools down. The partons can 
hadronize to mesons and baryons via inelastic interactions. When the chemical freeze-out 
point (temperature) is reached, the inelastic interactions will stop, and the relative parti- 
cle yields will not change. After that, the elastic interactions re-distribute the transverse 
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momenta among hadrons. The particle kinetic freeze-out point (temperature) will occur 
after the elastic interactions stop. 



1.2.1 Energy loss and jet quenching 

Particle yields will change due to energy loss in the strong interacting medium created 
in nuclear-nuclear collisions. These modifications of high px particle yields can be used as 
unique tools to study the medium properties. Experimentally, the widely used observable 
quantity of particle energy loss is called the nuclear modification factor Rab, which is 
defined as the ratio of the spectra in A + B collisions and those in p+p collisions, scaled 
by the number of binary nuclear- nuclear collisions |Adl02aj : 



^^^^^^^ - TABd^aPP/dprdv' ^^-^^ 

where Tab = {^bin) / o'inei accounts for the nuclear overlap geometry, averaged over the 
measured centrality class |Won94j . {Nbin) is the equivalent number of binary collisions 
calculated from the Glauber model. 

As presented in Fig. 11.51 (a), high px suppression of A^bin scaled hadron production 
in central Au+Au collisions relative to p+p collisions has been observed in RHIC ex- 
periments at y^s^=200 GeV. This is considered as evidence for the energy loss of the 
energetic partons, high px jets interacting with the hot and dense medium created in cen- 
tral Au+Au collisions at RHIC }Ada03b] . Due to the absence of nuclear effects such as 
shadowing |EKS981 IPFSOSl |Kop031 |Vogt04| , the Cronin effect |Accn41 1(:;atn41 IBPLFn4j . 



etc., hard processes are expected to scale with the number of binary collisions {Rab = 1)- 
Thus the nuclear suppression phenomenon was not seen in d+An collisions. The enhance- 
ment in the intermediate pt region for RdAu is due to the Cronin effect. The consistent 
picture of away-side jet quenching has also been observed in dihadron azimuthal angle 
correlations. Fig. 11.51 (b) shows the associated hadron {2<pj'<PT{trig) GeV/c) azimuthal 
distribution relative to a triggered hadron (4<j9r<6 GeV/c) after subtraction of the el- 
liptic flow and pedestal contributions. The enhanced correlation at near-side (A0 ~ 0), 
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which means the pair is from a single jet, was observed in p+p, d+Au and Au+Au 
coUisions. The pair from the back-to-back jet correlation at away-side (A0 ~ it) in cen- 
tral Au-|-Au collisions shows a dramatic suppression relative to those in p+p and d+Au 
collisions |Ada03bllW5l04j . 




Figure 1.5: Panel (a): Rab{pt) for minimum bias and central d+Au collisions and 
central Au-|-Au collisions. The bands show the normalization uncertainties. Panel (b): 
Two particle azimuthal distributions in p+p, d+An and central Au+Au collisions. 

All these high pj- suppression or away-side jet quenching observations indicate that 
light quarks strongly lose energy in the interacting dense medium. However, due to their 
large quark mass and small gluon radiative angle, heavy quarks are predicted to lose less 
energy compared to light quarks since pQCD energy-loss calculation assumes only gluon 
radiation |DK01l IDGW051 lArmOS] . The finite mass generalization of the small x (soft 
radiation x <^1) invariant DGLAP radiation spectrum is given by 

where uj{k) is the energy carried by radiated gluons in the medium with momentum fc, 
and Us is the strong coupling constant. The color charge factor Cr = = | with 

the number of color flavor A^^^ = 3 in this case. The gluon radiation angle is defined as 
6c = ^yrn^'+~x^M^ / (xE) . Due to large quark mass M, the gluon radiation is suppressed 
at smaller angles 6<6c- This effect is known as the "dead cone" phenomenon. 

From study on the non-abelian analog of the Ter-Mikayelian plasmon effect for gluons 



11 



5 10 15 20 25 

E (GeV) 



0.2 



5 10 15 20 

(GeV/c) 



Figure 1.6: Panel (a): The 1** order in opacity fractional energy loss for heavy quarks 
is shown as a function of their energy in a plasma characterized by as=0.3, Debye mass 
fi=0.5 GeV, and L = 5A = 5 fm. Panel (b): Heavy quark nuclear modifications as a 
function of pt for fixed L = 5 fm and dNg/dy = 1000. Short- (bottom) and long- 
(charm) dot-dashed curve include only radiative energy loss, while dashed (bottom) and 
solid (charm) curve include elastic energy loss as well. 

and the extended GLV energy loss method, the heavy quark energy loss to first order in 
the opacity, L/Xg {L is the path length, is the gluon mean free path), are numerically 
estimated as the solid curves (charm) and long-dashed curve (bottom) shown in Fig. 11.61 
(a), which is much less than light quarks (short-dashed curve). The shaded area indicates 
the so called "dead cone" effect. 

However, recent measurements of the pt distributions and nuclear modification factor 
for the non-photonic electrons (presumably from heavy quark decays) at high pt show 
similar suppression to those from light hadrons [STACRAAl IPHECRAA] , see Fig. 11.71 
This renews the interest in understanding the heavy quark energy loss mechanism in 
central Au+Au collisions at RHIC. From theoretical calculations, which take elastic (col- 
lisional) energy loss into account, the nuclear modifications of heavy quarks are expected 
to have larger suppression compared to those with only radiative energy loss |Wic05j . see 
Fig. 11.61 (b). Accounting for elastic energy loss seems to produce good agreement with 
experiments for heavy flavor energy loss, but this prescription fails in the case of light 
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Figure 1.7: Nuclear modification factors of non-pliotonic electron from PHENIX (open 
circle) and STAR (solid circle) are suppressed in similar magnitude with light flavor 

hadrons (grey band). 



flavors. These will be discussed in more detail in chapter 5. 



1.2.2 Collective motion and freeze-out 

The measured hadron spectra, especially in the soft sector at transverse momenta 
PT^l.d GeV/c, reflect the properties of the bulk of the matter produced in heavy-ion 
collisions after elastic interactions have stopped among the hadrons at kinetic freeze- 
out. At this stage the system is already relatively dilute and "cold". However, from the 
final state hadron spectra at kinetic freeze-out, one can obtain the information about 
the earlier hotter and denser stage. Since different hadrons have different production 
(hadronization) mechanisms, the characteristics of the different stages of the collision 
systems can be explored by analysis of the transverse momentum distributions for various 
hadron species. Fig. 11.81 (a) shows the measured tjit (= ^/pr + mass'^) spectra for light 
hadrons [n^K^p), A, S and multi-strange hadrons ($, Vt) in 200 GeV central Au+Au 
collisions |Ada04al IAda04bl IAda04cl ICas04j . and charmed hadron (D°) |YF06j in 200 
GeV minimum bias Au+Au collisions (will be discussed in Chapter 4). 
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Figure 1.8: Panel (a): The rriT spectra for light hadrons {7r,K,p), A, S and 
multi-strange hadrons ($, Q) in 200 GeV central Au+Au collisions, and charmed 
hadron {D^) in 200 GeV minimum bias Au+Au collisions are shown in symbols. The 
Blast Wave fit results are shown in curves. The BW fit for was done by combining 
spectrum and the measured charm decay leptons spectra below 2 GeV/c. The BW 
fits were done for vr", K~ , p simultaneously and for other particles separately. The 
arrows show the expected increasing freeze-out temperature and decreasing collective 
velocity from bottom to top. Panel (b): Blast-wave parameters Tfo vs. (/^r) contour 
plot from the simultaneous fits to light hadrons {7c,K,p) spectra in 5-10% and 0-5% 
Au+Au collisions and separate fits to multi-strange hadrons fi(sss) spectra in 

0-10% central Au+Au collisions. 
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In the Blast Wave thermal model |SSH93j . local thermal equilibrium is assumed and 
only particle mass and temperature of the system modify particle transverse momen- 
tum. After the collision, the bulk system expands and becomes more and more dilute 
and cold, while the particle collective velocity develops larger and larger. Extraction of 
model parameters characterize the random, generally interpreted as a kinetic freeze-out 
temperature Tfo and collective, radial flow velocity (/?), aspects. Fig. 11.81 (b) shows the 
fit parameters Tfo versus (jS) for different particles in central Au+Au collisions. Simul- 
taneous fit results for light hadrons (vr, K, p) give large flow velocity and small freeze- 
out temperature in central collisions, indicating that light hadrons kinetically freeze out 
late with strong collectivity when system evolutes a more rapid expansion after chemical 
freeze-out. Compared to light hadrons, the larger temperature and smaller flow velocity 
obtained from the flts indicates that multi-strange hadrons freeze out at earlier stage dur- 
ing the evolution of the bulk system. Since multi-strange hadrons are expected to have 
smaller hadronic scattering cross sections and their transverse momentum distributions 
will not change signiflcantly after chemical freeze-out, the kinetic freeze-out temperature 
from the flt to those particles is consistent with the chemical freeze-out temperature Tch, 
albeit with still large uncertainties. The chemical freeze-out temperature Tch is close to 
the critical temperature Tc, indicating the temperature of the system created in the col- 
lisions is greater than Tc and hence the phase transition may take place and QGP may 
form at RHIC energy of ys^=200 GeV. 




Due to relatively heavy quark mass and much smaller hadronic scattering cross section, 
heavy flavor hadrons are expected to freeze-out early and difficultly participate in collec- 
tive motion. Thus the larger freeze-out temperature and smaller fiow velocity are expected 
for heavy fiavor hadrons. But currently the experimental statistics is not good enough 
to distinguish the freeze-out properties between D-meson and multi-strange hadrons, this 
will be discussed in chapter 5. 
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1.2.3 anisotropic flow 



Recent heavy-ion collisions experiments at RHIC have successfully generated a lot of 
exciting results. Besides strong suppression of particles with high momentum created in 
the collisions - jet quenching, the large azimuthal anisotropy - elliptic flow is one of the 
most remarkable discoveries. All observations indicate that the extremely opaque quark- 
gluon matter formed in the collisions exhibits highly collective, near-hydrodynamic behav- 
ior |KH03j . In non-central heavy- ion collisions, the initial coordinate spacial anisotropy, 
which is formed by the geometrical overlap of the two nuclei, will produce anisotropic 
pressure gradients in the high dense matter which then subsequently produce a final-state 
momentum space anisotropy of the produced particles via rescattering. Because of the 
self-quenching of coordinate space anisotropy, the measurement of the particle azimuthal 
anisotropy distributions with respect to (w.r.t) the reaction plane can reveal the informa- 
tion on the dynamics at the early stage of the collision. 

The initial space anisotropy in the overlapping region (x,y), which is the plane trans- 
verse to the beam axis z, is characterized by the eccentricity: 

(y^ + x^) ^ ^ 

The final particle spectrum in momentum space can be expanded into a Fourier ex- 
pression in terms of particle azimuthal distributions w.r.t the reaction plane '^r'- 

The Fourier coefficients represent anisotropy parameters and can be extracted as 

t;„ = (cos K0 (1.12) 

where the first and second coefficients fi, V2 are called the directed and elliptic flow. 
The identifled particle elliptic flow f 2 shows strong transverse momentum dependence. 
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In the low pt region, particle collective motion shows hydrodynamical behavior |KH03t 
lAdaOSj . As is well known, hydrodynamic models, which assume ideal relativistic fluid 



flow and negligible relaxation time compared to the time scale of the equilibrated system, 
successfully reproduce the experimental data and reasonably describe the mass ordering 
for V2 at low pt- Thus the particle V2 follows ttit — m scaling |Yan07] . as shown in 
Fig. 11.91 (a). Here ttit = \/p\ ^ irfi and m is the particle mass. On a m-r — m scale, 
at intermediate pt (1-5 GeV/c < pt < 4: GeV/c), the particle V2 is split into meson and 
baryon groupings, and this phenomenon may be described by coalescence/recombination 
models |LKn2l[MVn3] . In these models, mesons and baryons are hadronized by coalescing 
two and three co-moving quarks, respectively. The simple picture is that the hadron V2 
can be described by its constituent quark t>2 as: 

V2{Pt) — ngV2{pT / n-q) , nq{meson) = 2, nqiparyon) = 3, (1-13) 

where ng is the number of constituent quarks. Thus, V2 may follow the Number- of- Quark 
(NQ) scaling. This scaling behavior is observed in Au+Au collisions at RHIC. Fig. 11.91 
(b) shows the scaling of identified particle f 2 divided by function of (mj- — m)/nq. 

These exciting results are considered as indicators for the observation of the partonic level 
collectivity (deconfinement) in heavy-ion collisions at RHIC energies. 



1.2.4 Charm production at RHIC 

Charm quarks are a unique tool to probe the partonic matter created in relativistic 
heavy-ion collisions at RHIC energies. Because of their large quark mass 1.3 GeV/c^), 
charm quarks are predicted to lose less energy compared to light quarks when pQCD 
energy-loss calculation assumes only gluon radiation |DK01l IDGW051 lArmOS] . Recent 
measurements of the pt distributions and nuclear modification factor for the non-photonic 
electrons (presumably from heavy quark decays) at high p^ show similar suppression to 
those from light hadrons |STACRAA[ IPHECRAAj . This renews the interest in under- 
standing the picture of the observed jet-quenching phenomena at the partonic stage in 
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Figure 1.9: Panel (a): Identified particle V2 as a function of — m in 0-80% Au+Au 
collisions at ys^=200 GeV. Open circles, open squares, solid triangles and solid circles 

represent -ft^g , A + A, S + S and + f2, respectively. Hydrodynamic calculations are 
shown as curves for comparison. Panel (b): Identified particle scaled by the number 
of constituent quarks (n^) versus (m^ — m)/nq. 

Au+Au collisions at RHIC Furthermore, theoretical calculations have shown that inter- 
actions between the surrounding partons in the medium and charm quarks could change 
the measurable kinematics [MT051 IHGRO^ IDDZ06] . and could boost the radial and ellip- 
tic flow resulting in a different charm px spectrum shape. Due to its heavy mass, a charm 
quark can acquire flow from the interactions with the constituents of a dense medium in 
analog to Brownian motion. Its decoupling from medium also depends on the cross-section 
of the interaction and medium density. Therefore, the measurement of charm flow and 
freeze-out properties is vital to test light flavor thermalization and the partonic density 
in the early stage of heavy ion collisions pOlOSl ILM031 IGKR041 IMT051 IBKG03j . The 
charmed hadron kinetic freeze-out temperature Tfo and the radial flow velocity 13 can be 
derived from blast-wave model |SSH93] fits to data. 

These postulations and predictions are based on the assumptions that charm quarks 
are produced only at early stages via initial gluon fusion and that their production cross- 
section can be calculated by perturbative QCD [ Lin95t ICNVOSj . Study of the binary 
collision (Abin) scaling properties for the total charm cross-section from d+Au to Au+Au 
collisions can test these assumptions and determine if charm quarks are indeed good 
probes with well-defined initial and final states. Charm total cross-section measurement 
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is also essential for model calculations, which tries to explain observed similar suppression 
pattern of J/^ at RHIC and SPS [PbmOOl IABR031 IGROTI IKos02[ ITSROli IGKL02j . It 
has been shown that low pt muons and intermediate Pt electrons from charmed hadron 
semileptonic decays can stringently constrain the charm total cross-section and are sensi- 
tive to the radial flow and the freeze-out condition |LZZX06] . 

1.2.5 Bottom contribution in non-photonic electron measurements 

Due to the absence of the measurement of B-mesons and precise measurement of 
D-mesons, it is difficult to separate bottom and charm contributions experimentally in 
current non-photonic electron measurements for both spectra and elliptic flow V2- As 
discussed previously, the suppression behavior of heavy quarks is quite different from 
light quarks due to the "dead cone" effect |DK01j . and this is especially true for the 
bottom quark. Even when the elastic energy loss is included, the bottom quark still loses 
much less energy. The bottom contribution may reduce the energy loss of non-photonic 
electrons from heavy flavor decays. But as shown in Fig. 11.101 (a), the suppression of 
the non-photonic electron Raa is as large as light hadrons. Both the theoretical result 
with charm energy loss only and the theoretical calculations with charm+bottom energy 
loss by assuming large q or counting elastic energy loss can describe the data within 
errors |DGV06l IWicOSl [RapOGj IVit06| lArmOGj . 



Recently, PHENIX has measured the non-photonic electron V2 |SS07j . shown in Fig. II. 10] 
(b). The observed large elliptic flow of the non-photonic electron may indicate strong cou- 
pling of heavy quarks with medium. There are many theoretical calculations for the non- 
photonic electron such as charm thermal+fiow model (solid curve) |GKR04j . A multi- 
phase transition (AMPT) model assume cross section (jp=10 mb (dashed curve) |ZCK05j . 
resonance states of D-/B- mesons (band) jHGROGj . etc. The comparison with theories also 
showes that both the model results with charm only and the results with charm+bottom 
have good agreement with data within errors. 

Thus, the puzzle of the bottom contributions in non-photonic electron spectra and V2 
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Figure 1.10: Panel (a): Non- photonic electron Raa measured from STAR (circle) 
compared to theories. Panel (b): Non-photonic electron V2 measured from PHENIX 

(circle) compared to theories. 



still remains. We will have a detail discussion about this issue in chapter 5. 
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CHAPTER 2 



Experimental Set-up 

2.1 The RHIC accelerator 

The Relativistic Heavy Ion Collider (RHIC) located at Brookhaven National Labora- 
tory (BNL) is built to accelerate and collide heavy ions and polarized protons with high 
luminosity. The RHIC accelerator facility consists of two concentric storage rings, called 
blue and yellow rings, sharing a common horizontal plane in the tunnel. The super- 
conducting magnets in each ring with a circumference of 3.8 km are designed to bend 
and focus the ion beams. The counter-rotating beams can collide with one another at 
six location along their 3.8 km circumference. The top center-of-mass collision energy 
for heavy ion beams is 200 GeV per nucleon pair. The operational momentum increases 
with the charge-to-mass ratio, resulting in the top energy of 125 GeV/u for lighter ion 
beams and up to 250 GeV/u for polarized proton beams. The average luminosity for 
gold-on-gold collisions at y/s^= 200 GeV is 8 x 10^^cm~^s~^ without electron cooling 
and 7 x 10^'^ cm~^s^^ with electron cooling. 

Fig. 12.11 shows the RHIC accelerator complex including the accessorial accelerators 
used to bring the ions beams up to RHIC injection energy, and strip electrons from the 
atoms. Negatively charged gold ions are partially stripped of their electrons and then 
accelerated to 15 MeV/u in the Tandem Van de Graaff facility. After a charge selection 
by bending magnets, beams of gold ions are transferred to the Booster Synchrotron and 
accelerated to 95 MeV/u through the Tandem-to-Booster line. Then the gold ions are 
injected into the Alternating Gradient Synchrotron (AGS) and accelerated to 10.8 GeV/u. 
The beams are transferred to RHIC through the AGS-to-RHIC Beam Transfer Line. 
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Finally, beams are injected to RHIC and fully accelerated to the top collision energy 100 
GeV/nucleon. For p+p collisions, proton beams are injected from the 200 MeV Linac 
into the booster, followed by acceleration in the AGS and injection into RHIC. 




Figure 2.1: Schematic of the RHIC complex. RHIC's two 3.8-kilometer rings collide 
relativistic heavy ions and polarized protons at six intersection points. 

There are six interaction points along the rings, and 4 of them are equipped with 
detectors. They are two large experiments STAR (6 o'clock), PHENIX (8 o'clock) and 
two small ones PHOBOS (10 o'clock) and BRAHMS (2 o'clock), respectively. 

2.2 The STAR detector 

To investigate the behavior of strongly interacting matter at high energy density and 
to search for signatures of the new matter form of QGP, The Solenoidal Tracker at RHIC 
(STAR) was specially constructed for measurements of hadron production over a large 
solid angle with high precision tracking, high quality momentum, and good particle iden- 
tification capability. It has an azimuthal symmetric acceptance and covers large range 
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around mid-rapidity. STAR consists of several subsystems and a main tracker - the Time 
Projection Chamber (TPC) located in a homogenous solenoidal analyzing magnet. 




Figure 2.2: Perspective view of STAR detector, with a cutaway for viewing inner 

detector systems. 

The layout of the STAR experiment is shown in Fig. 12. 2[ Fig. 12.31 shows the cutaway 
view of the STAR detector. The STAR magnet is cylindrical in design with a length 
of 6.85 m and has inner and outer diameters of 5.27 m and 7.32 m, respectively. A 
uniform magnetic field of 0.5 T (full field or reversed full field) or 0.25 T (half field) 
is provided for the tracking of charged particles. The main tracker - TPC for charged 
particle tracking and identification is 4 meters long and located at a radial distance from 
50 to 200 cm from the beam axis including 45 layers. It covers a pseudo-rapidity range 
\ri\ < 1.5 and in full azimuth (27r). There are inner detectors Silicon Vertex Tracker 
(SVT) and Silicon Strip Detector (SSD) close to the interaction vertex, which provides 
additional high precision space points on track so that it improves the position resolution 
and allows us to reconstruct the secondary vertex of weak decay particles. Both the TPC 
and SVT contribute to particle identification using ionization energy loss (dE/dx), with 
an anticipated combined energy loss resolution of a = 7%. The momentum resolution 
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reaches a value of 5p/p — 0.02 for a majority of the tracks in the TPC. The momentum 
resolution improves as the number of hit points (nhitpts) along the track increases and 
as the particles momentum decreases, as expected. 
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Figure 2.3: Cutaway view of the STAR detector. It includes the partial installed 
ElectroMagnetic Calorimeter (EMC) and two prototypes Time-of-Flight (TOP) 

detectors. 



In the extended large rapidity region, there are two Forward TPC (FTPC) detectors 
installed covering 2.8 < I77I < 3.8 with complete azimuth to track particles at forward 
and backward rapidity. One prototype tray of Time- Of -Flight (TOF) detector using 
scintillator materials (TOFp) was installed since Run II and another prototype tray of 
TOF detector using Multigap Resistive Plate Chamber (MRPC) technology (TOPr) was 
installed since Run III, covering — 1 < jjyl < and ^ of full azimuth. In future run, 
TOFp will be removed and a full barrel TOP detector based on MRPC technique will 
be installed, which is expected to extend the PID capability of STAR greatly. Part 
of barrel ElectronMagnetic Calorimeter (BEMC) and endcap EMC (EEMC) were also 
installed since Run II. They are used to measure the electromagnetic probes - electrons 
and photons. 

The fast detectors are designed for providing inputs of the trigger system, such as 
Zero-degree Calorimeter (ZDC) , Central Trigger Barrel {CTB) and Beam-Beam Counters 
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(BBC). Two ZDCs locates on each side ~ 18 m away from the colhsion points. Each 
is centered at 0° and covers ~ 2.5 mrad. In the colhsion, usually just the overlap parts 
of the two colliding ions interact with each other. The left fragments are continuously 
going forward. The charged fragments are bent by the dipole magnets and the outgoing 
neutrons can be detected by the ZDCs. The ZDC signals are used for monitoring the 
heavy ion beam luminosity and for the experiments triggers. A minimum bias trigger was 
obtained by selecting events with a pulse height larger than that of one neutron in each 
of the forward ZDCs, which corresponds to 95 percent of the geometrical cross section. 
The CTB is a collection of scintillating tiles surrounding the outer cylinder of the TPC. 
The CTB will be mostly used to select central triggered events in heavy ion collisions by 
measuring the occupancy of those CTB slats. The BBC subsystem covers 3.3 < \ri\ < 5.0, 
measuring the "beam-jets" at high rapidity from Non- Singly Diff'ractive (NSD) inelastic 
p+p interactions. It consists of two disk shaped scintillating detectors, with one placed 
at each endcap of the TPC (3.5 m from TPC center). Each BBC disk is composed of 
scintillating tiles that are arranged in a hexagonal closest packing. The j>+j>NSD trigger 
sums the output of all tiles on each BBC and requires a coincidence of both BBC's firing 
above noise threshold within a time window. 



2.3 The Time Projection Chamber - TPC 

TPC is the primary tracking device of the STAR detector |And03] . The TPC is 
designed to record the tracks of particles, provide the information of their momenta, and 
identify the particles by measuring their ionization energy loss [dE/dx). Charged particles 
are identified over a large momentum range from 0.15 GeV/c to ~ 30 GeV/c. For half 
field, the lower limit ~ 0.075 GeV/c can be achieved. The TPC consists of a 4.2 m long 
cylinder with 4.0 m in diameter. The cylinder is concentric with the beam pipe, and the 
inner and outer radii of the active volume are 0.5 m and 2.0 m, respectively. The TPC 
covers the full region of azimuth (0 < < 27r) and covers the pseudo-rapidity range of 
\ri\ < 2 for inner radius and \rj\ < 1 for outer radius. Fig. 12.41 shows a cutaway view of the 
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TPC structure. 




Eld 



Figure 2.4: Cutaway view of the TPC detector at STAR. 



The TPC is divided into two parts by the central membrane. It is typically held at 
28 kV high voltage. A chain of 183 resistors and equipotential rings along the inner and 
outer field cage create a uniform drift filed (~ 135 V/cm) from the central membrane to 
the ground planes where anode wires and pad planes arc organized into 12 sectors for 
each sub-volume of the TPC. The working gas of the TPC is PIO, the mixture of 90% Ar 
and 10% CH4, regulated at 2 mbar above the atmospheric pressure. The electron drift 
velocity in PIO is relatively fast, ~ 5.45 cm/fis at 130 V/cm drift field. The gas mixture 
must satisfy multiple requirements and the gas gains arc ~ 3770 and ~ 1230 for the 
inner and outer sectors working at normal anode voltages (1170 V for inner and 1390 
V for outer), respectively. The readout system is based on the Multi-Wire Proportional 
Chamber (MWPC) with readout pads. Each readout pad is also divided into inner and 
outer sub-sectors, while the inner sub-sector is designed to handle high track density 
near collision vertex. 136,608 readout pads provide {x, y) coordinate information, while 
z coordinate is provided by 512 time buckets and the drift velocity. Typical resolution is 
~ 0.5 — 1.0 mm. 
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When charged particles traverse the TPC, they hberate the electrons from the TPC 
gas due to the dE/dx. These electrons are drifted towards the end cap planes of the 
TPC and collected by a readout pad. The signal is amplified and integrated by a circuit 
containing a pre-amplifier and a shaper. Then it is digitalized and then transmitted over 
a set of optical fibers to STAR Data AcQuisition system (DAQ). 

At the DAQ stage, raw events containing millions of ADC values and TDC values 
were recorded. Raw data were then reconstructed into hits, tracks, vertices, the collision 
position through the reconstruction chain of TPC by Kalman method. The TPC recon- 
struction process begins by the 3D coordinate space points finding. This step results in 
a collection of points reported in global Cartesian coordinates. The Timing Projection 
chamber Tracker (TPT) algorithm is then used to reconstruct tracks by helical trajectory 
fit. The resulted track collection from the TPC is combined with any other available 
tracking detector reconstruction results and then refit by application of a Kalman filter 
routine — a complete and robust statistical treatment. The primary collision vertex is 
then reconstructed from these global tracks and a refit on these tracks with the distance 
of closest approach (dca) less the 3 cm is preformed by a constrained Kalman fit that 
forces the track to originate from the primary vertex. As expected, the vertex resolution 
decreases as the square root of the number of tracks used in the calculation. The primary 
vertex resolution is ~ 350 fim with a track multiplicity above 1000. The reconstruction 
efficiency including the detector acceptance for primary tracks depends on the particle 
type, track quality cuts, pt, track multiplicity etc. The typical value for the primary pions 
with Nfit > 25 and \ri\ < 0.7, dca < 3.0 cm is approximate constant at pt > 0.4 GeV/c: 
>~ 90% for Au+Au peripheral collisions and ~ 80% for central collisions, respectively. 

The TPC provide the track momentum and the dE/dx information for charged parti- 
cles identification. For a particle with charge z (in units of e) and speed [3 = v/c passing 
through The mean rate of dE/dx is given by the Bethe-Bloch equation 12.11 [PDG04j : 



dE 
dx 




max 



2 



(2.1) 
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The meaning of each symbol can be referred to |PDG04] . Different types of particles 
(different rest masses) with the same momentum have different kinematic variables (3 
(7), which may result in distinguishable dE/dx. The typical resolution of dE/dx in 
Au+Au collisions is ~ 8%, which makes the n/K separation up to p ~ 0.7 GeV/c and 
proton/meson separation up to p ~ 1.1 GeV/c. 

Combined the TPC with other detectors, such as the TOF detector and the BEMC, 
etc., the capability of particle identification can be greatly improved, particle separation 
can be extended to higher region. 

2.4 The prototype TOF detector based on MRPC technology 

A full barrel TOF detector, based on the MRPC technology, has been proposed in the 
STAR detector upgrade programs. There will be 120 trays with 60 on east side and 60 
on west side. For each tray, there will be 33 MRPCs. The TOFp detector (a prototype 
based on scintillator technology) was installed since Run II |Llo04] . It replaced one of 
CTB trays, covering —1 < < 0, and 7r/60 in azimuth. It contains 41 scintillator slats 
with the signal read out by Photo Multiplier Tubes (PMTs). The resolution of TOFp is 
~ 85 ps in Au+Au collisions. However, due to the significant higher cost by the PMTs, 
this design will not be used in the full TOF upgrade. 

On the other hand, in Run III and Run IV, new prototypes of TOF detector based 
on MRPC (TOFr) were installed also covering — 1 < 77 < and 7r/60 in azimuth. In Run 
III, 28 MRPC modules were installed in the tray and 12 of them were equipped with 
electronics, corresponding to ~ 0.3% of the TPC acceptance |Ada03cj . In Run IV, a new 
tray with 24 modules were installed at the same place as Run III. But only 12 modules 
were equipped with valid electronics, which means the acceptance in Run IV was roughly 
similar to that in Run III. 

The trigger system of the TOF detector is the two pVPDs, each staying 5.4 m away 
from the TPC center along the beam line |Llo04j . They provide a starting timing infor- 
mation for TOF detectors. Each pVPD consists of three detecting element tubes covering 
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~ 19% of the total solid angle in 4.43 < \t]\ < 4.94. Due to different multiplicities, the 
effective timing resolution of total starting time is 25 ps, 85 ps and 140 ps for 200 GeV 
Au+Au, d+Au and p+p collisions, respectively. 
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Figure 2.5: Two-side view of a MRPC module |Sha02] . 

Resistive Plate Chambers (RPCs) were developed in 1980s. MRPC technology was 
first developed by the CERN ALICE group less than 10 years ago. Fig. 12.51 shows the two 
side views (long edge view on top and short edge view on bottom) of an MRPC module 
appropriate for STAR |Sha02j . An MRPC basically consists a stack of resistive plates, 
spaced one from the other with equal sized spacers (such as fish line) creating a series 
of gas gaps. It works in avalanch mode. Electrodes are connected to the outer surfaces 
of the stack of resistive plates while all the internal plates are left electrically floating. 
Initially the voltage on these internal plates is given by electrostatics, but they are kept 
at the correct voltage due to the flow of electrons and ions created in the avalanches. 
There are six read-out strips on each module in this design. The first beam test for 6-gap 
MRPCs at CERN PS-TIO facility with piab = 7 GeV/c pions beam resulted in a ~ 65 ps 
timing resolution with more than 95% detecting efficiency and the module is capable of 
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working at high event rate (500 Hz/cm^) |Sha02j . These modules were then assembled 
in a prototype TOF tray and tested in the AGS radiation area. Similar resolution was 
obtained. In RHIC Run III and Run IV, the MRPC modules in TOFr trays installed in 
the STAR detector were applied on the high voltage of 14 kV and with the working gas of 
95% freon and 5% isobutane. The charged particle detecting efficiency is > 95% at high 
voltage plateau. 

TOF system calibrations include the start time calibration from pVPDs and TOFr/TOFp 
flight time calibration. The main sources need to be considered are global time offset due 
to different electronics delays, the correlation between the amplitude and the timing 
signals, the correlation between the hit position and the timing signals etc. Detailed 
calibrations on TOF systems can be found in |Rua04bt IAda03c] (TOFr) and |Llo04] 
(TOFp). 
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CHAPTER 3 



Non-photonic electron measurement in 200 GeV 

Au+Au collisions 

Due to the large combinatorial background, direct reconstruction of heavy-flavor 
mesons via hadronic decay channels is diflicult in current high energy nuclear collisions 
at RHIC jSTAcOSj . The indirect measurement of single leptons through heavy-flavor 
semileptonic decays has been used as an efficient way to study heavy flavor production. 
The transverse momentum (px) distributions of single electron and the extracted total 
charm cross-section have been measured in 200 GeV d+An collisions |STAc05[ IXinOSj . 
Study of the binary collisions (A/'bin) scaling properties for the charm total cross-section 
among d+Au to Au+Au collisions can test if heavy-flavor quarks are produced exclusively 
at initial impact. Modification of the single lepton production in nuclear collisions can re- 
veal the heavy-flavor energy-loss. Thermal parameters extracted from the px distributions 
of single electron are useful for us to understand the heavy-flavor hadrons freeze-out and 
flow properties in medium in nuclear-nuclear collisions. Therefore, it is very important to 
measure single electron px distributions in Au+Au collisions. 

In this chapter, analysis details of single electron pt distributions in 200 GeV Au+Au 
collisions will be presented. In addition, due to large photonic background, it is difficult to 
measure single electron V2 in current experimental environment, but a method to extract 
single electron t>2 will be proposed. 
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3.1 Non-photonic electron transverse momentum distributions 

3.1.1 Data sets and cuts 

In RHIC Run IV (year 2004), STAR experiment collected abundant triggered mini- 
mum bias (minbias) events and central events from y/sNN=200 GeV Au+Au collisions. A 
0-80% minbias Au+Au collision centrality is defined by the charged particle reference mul- 
tiplicity (refmult) measured at mid-rapidity (|?7| < 0.5) requiring number of TPC fit points 
(nFitPts) > 9 and the distance-of-closest-approach (dca) < 3 cm. The sub-centralities 
for the data analysis are selected by applying the refmult cuts from the calculation of 
Glauber model [GlauOlj . which gives the average percentage of the number of events in 
each centrality bin, see Fig. 13. 1[ For the centrality dependence study of the non-photonic 
electrons, the minbias event sample was subdivided into three centrality bins: 0-20%, 
20-40% and 40-80%. Table \3A\ lists the data sets under the trigger selections used in this 
analysis from Run IV 200 GeV Au+Au collisions. The right column lists two numbers 
for the size of event samples after the location of collision vertex along beam axis (V^) 
cuts, one is the number of TPC triggered events used for the photonic background elec- 
tron analysis, the other is the number of TOF triggered events used for inclusive electron 
identification. 



Table 3.1: Data sets from Run IV used in this analysis 



Trigger 


Centrality(refmult cut) 


Vertex Z cut (cm) 


Events Size (TPC, TOF) 


minbias 


0-20% [319,1000) 




Vz 


< 30 


3.28M, 1.93M 


20-40% [150,319) 




Vz 


< 30 


3.36M, 1.97M 


40-80% [14,150) 




Vz 


< 30 


6.55M, 3.85M 


0-80% [14,1000) 




Vz 


< 30 


13.2M, 7.8M 


central 


0-12% 




Vz 


< 50 


21. 7M, 15.5M 



The online centrality selection in Au+Au minbias collisions was based on the energy 
deposited in the two Zero-Degree Calorimeters (ZDC) |STA03l IAda05j . Since the 5- 
10% refmult cuts in central triggers will give some combination of ZDC- and TPC-based 
centrality, which is not 5-10% actually. But in the analysis, due to the limited statistics, 
all the central triggers are accepted, the centrality is 0-12%, the bias from realistic 0-10% 
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is about 5%. The TOF trigger set up in Run IV was to select events with a vahd pVPD 
coincidence. The trigger efficiency for the TOF events is ~60%. The tracks matched to 
the TOF should have good ADC signal (> 30) and at least one TOF hit. In Run IV, one 
TOFr tray, which is based on MRPC technology, and one TOFp tray, which uses RFC 
modules, were installed outside the TPC and calorimeters. The TOF detector covers 
1/120 full barrel and negative pseudorapidity (— 1 < ?7 < 0). 

10^ 

w 
c 

o 
O 

10^ 



1 

200 400 600 800 

Multiplicity 

Figure 3.1: Panel (a) shows the multiplicity for minbias triggered events, which follows 
Glauber distribution. Panel (b) shows the centralities by cutting the multiplicity 

calculated from Glauber model. 

3.1.2 Electron identification 

The charged particle is usually identified by the ionization energy loss (dE/dx) mea- 
sured in the STAR main detector TPC |And03] . With TPC dE/dx only, electron can be 
identified only up to 0.8 GeV/c |Ada04dl IJoh02j due to large hadron contaminations. 
But the very small electron mass gives it the ability to be separated from other hadrons 
by measuring the flight timing information. The TOF detector was designed to measure 
the velocity of electrons and hadrons [TofOSl ISTAcOSl IShaOGl IXinOSj . 

By using a combination of velocity measured from the TOF and dE/dx measured 
in the TPC, inclusive electrons are identified up to 4 GeV/c in minbias Au+Au 
collisions and up to 5 GeV/c in central Au+Au colhsions. Panel (a) of Fig. 13.21 
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shows the 2-D scattering plot of dE/dx measured from TPC vs. momentum (p) for the 
charged particles with good TOF hits matched in Au+Au collisions. Panel (c) shows the 
l/[3 vs. p. Due to the different mass, charged particles can be separated by measuring 
their velocity. Electron band is around unity and merged with pion band. Panel (b) 
shows the dE/dx vs. p after the particle velocity cut — 1| < 0.03), after most of the 
charged hadrons rejected, the pure dE/dx bands for electron and pion are left. And they 
are separated well. 




12 3 

Momentum p (GeV/c) 



Figure 3.2: Panel (a): Ionization energy loss measured by TPC vs. momentum. Panel 
(b): Ionization energy loss measured vs. momentum after TOF velocity cut. Panel (c): 
Particle velocity measured by TOF vs. momentum. 

Table 13.21 lists the cuts for inclusive electron selection. 

The dE/dx distributions are not exactly gaussian. The logio{{dE/dx)/{dE/dx^-^f^^^i)) 
(abbreviated as dE/dx), which follows gaussian distribution, was defined as the logarith- 
mic ratio of measured dE/dx and predicted dE/dx from Bichsel model [BicOGj . After the 
velocity cut, the 2-dimension histogram of dE/dx as a function of pr was projected into 
several pt bins. The 2-Gaussian function cannot describe the shoulder region between 
hadron peak and electron peak well at lower px, a function of exponential+gaussian was 
performed instead. At GeV/c, the other hadrons, such as protons, contaminate 
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Table 3.2: Cuts for inclusive electron selection 



Primary track? 


yes 






Global dca 


< 3.0 cm 


nFitPts 


> 25 


nFitdedxPts 


> 15 


nFitPts/nMax 


> 0.52 




(-0.03, 0.03) 


pseudorapidity 


(-1.0, 0) 


TOFr ADC 


> 30 



the low dE/dx window and statistics cannot enable us to distinguish the difference of 
these two fits, so 3-Gaussian fit was used in this pr region. Then the inclusive elec- 
tron raw yields were obtained from these fits in each px bin. Fig. 13.31 shows the dE/dx 
distributions and fitting results in several typical pt bins after TOF velocity cut. The 
integral yields by applying a cut on 3(7 > dE/dx > dE/dx^^^^ are used to estimate the 
hadron contamination compared to the fits, dE/dx^^^^ is the central value of the dE/dx 
distribution. During the energy deposition in the TPC, the merged tracks are measured 
as larger energy loss, which dominate the tails shown on the right of electron gaussian 
peak, see Fig. 13. 3[ It is more obvious in central collisions than peripheral. The higher 
the multiplicity is, the larger the possibility of the track mergence. There are two small 
gaussian contributions to the tail: one is due to two merged pions, the other is considered 
as a mergence of one electron and one pion. The contamination of the tail to the electron 
is quite small, less than percent, even it looks obvious in the log scale. The tails fade at 
higher pr- 

Fig. 13.41 show the difference of the fit to the tails by fixing width and opening width. 
The yields of the inclusive electron change less than 1%. 

Since dE/dx projected in a momentum bin is more like gaussian distribution than 
that projected in a pt bin. The naipx) is defined as 

na{pT) ~ A X {dE/dx - dE / dx j^.^^^J / {dE / dx ^^^,^^^1 x VndEdxEitPts) (3.1) 
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Figure 3.3: dE/dx projections in several pt bins after the TOF velocity cut. The 
electron raw yields were extracted from the fit. Exponential (red curve) + gaussian 
(blue curve) function is used at low p^- Multi-gaussian function is performed at higher 
Pt- The tails due to merged tracks are described by gaussian functions, they fade at 
higher p^- The dE/dx distributions in minbias collisions (first 4 panels) and central 

collisions (last 2 panels) are similar. 
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Figure 3.4: Left panel shows the gaussian fit to the tail with fixed width. Right panel 
shows the fit with open width. The difference of electron raw yields is < 1%. 

y4 is a constant, which only depends on the detector environment, na is gaussian 
distribution in each bin. ndEdxFitPts is the number of fit points used for calculating 
the dE/dx of a charged particle during the tracking in the TPC. A 2-Gaussian function 
was used to fit the na^ shown as Fig. 13. 5[ The difference between fit to ncTg and fit to 
log dE/dx is less than 5% point to point at lower pt- At higher px the difference is even 
smaller, due to the smaller difference of momentum and pt- 
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Figure 3.5: The difference between fit to na^ and log dE/dx is less than 5% point to 
point at lower p^- This was used to check the difference of the electron raw yields 

derived from dE/dx vs. p and pt- 
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3.1.3 Acceptance and efficiencies 

During the detector tracking and tlie extraction of ttie inclusive electron raw yields, a 
part of them were lost due to the detector acceptance, the reconstruction and cut efficiency 
etc.. The final yields of inclusive electron need to be corrected by the acceptance and 
efficiencies: 

• TPC tracking efficiency — track reconstruction efficiency, TPC acceptance. 

• TOF matching efficiency — TPC track with good TOP hits matched, TOP accep- 
tance. 

• number of dE/dx fit points cut efficiency. 

Due to the lower statistics of electrons, the TPC tracking efficiency was studied from 
pion embedding for each centrality in Au+Au collisions, shown in left panel of Pig. 13.61 
The different vertex Z cuts and eta cuts do not effect the tracking efficiency much. The 
TOP matching efficiencies used for electrons were determined from real data as the number 
of pion tracks matched to TOP divided by the number of pion tracks reconstructed from 
TPC. The centrality dependence of the matching efficiency is shown in right panel of 
Pig. 13.61 But different from electrons, ~10% of pions could be scattered or decay, especially 
at low which was corrected in the total efficiencies. 

The TPC tracking efficiency does not include the ndEdxPitPts cut efficiency, since it is 
technically difficult to study dE/dx in embedding data. So the ndEdxPitPts cut efficiency 
was determined from real data as the number of tracks by cutting on ndEdxPitPtsgfeglS 
divided by the number of total tracks, shown in Pig. 13.71 

3.1.4 Photonic background 

In this analysis, the dominant sources of photonic electron background are considered 
as gamma conversions and scalar meson Dalitz decay: 

• 7 — i> e^e~ photon conversions in the material in STAR detector. 



38 



1 



>- 0.8 
c 

I 0.6 

Q 

J 0.4 

o 
n3 



0.2 




IVzl<20 cm, hikO.7 
IVzl<30 cm, hil<l 



minbias Au+Au 200 GeV 



40-80% 
20-40% 
0-80% 
0-20% 
0-10% 



2 4 

(GeV/c) 



o 

03 
O 



o 



0.5 




minbias Au+Au 200 GeV 



0-80% 
0-20% 
20-40% 
40-80% 



(GeV/c) 



Figure 3.6: Left panel: Centrality dependence of TPC tracking efficiency from 
embedding data. Right panel: Centrality dependence of TOF matching efficiency from 

real data. 
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Figure 3.7: Efficiency of cutting on ndEdxFitPts > 15 as a function of pt- 
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• 7r°^7e+e- (1.198 ± 0.032)%. 

• r]^ 7e+e- (0.60 ± 0.08)%. 

The huge photonic background electrons come from conversions, especially at low pt{< 
1 GeV/c), due to the amount of materials in STAR detector. Fig. l3.8l shows the conversion 
point and conversion radius distributions from the simulation data with HIJING generator 
+ STAR GEANT detector configuration in 200 GeV Au+Au collisions. It shows most of 
converted electrons are from the huge material of 3-layer SVT {R ~ 7, 11, 16cm) and its 
supporting structures. 



^ ^ ^ x10^ 




Conversion Point X (cm) Conversion Radius (cm) 

Figure 3.8: Left panel: The 2-D scattering plot for conversion point, where the 
photonic electron come from. It reflects the STAR material structure. Right panel: The 

conversion radius distributions. 



Photonic electrons dominate the total yield, especially at low pt- At pt> 1 GeV/c, 
The electrons from heavy flavor hadron decays, which is considered as the signal, may 
become visible due to the enhancement of the total spectrum w.r.t background spec- 
trum. From the previous measurements, background from photon conversion etc. was 
reconstructed experimentally using topological method |Ada04d[ IJoh02j . The e'^e" from 
photonic background have the characteristic feature of low invariant mass. The e~^e~ 
invariant mass were reconstructed from a tagged electron (positron) from TOF at low pt 
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(we use TPC dE/dx to identify electron at high pt-, due to the poor statistics from TOF 
events at high px) combined with every other global positron (electron) candidate, whose 
helixes have a dca (rfcag+e-) to the tagged tracks' less than 1 cm, to find the other partner 
track reconstructed in the TPC. Because the TPC acceptance is large enough, the pair 
reconstruction efficiency is reasonable. Two methods to reconstruct the invariant mass 
for e^e~ pairs were used in this analysis: One is to reconstruct the invariant mass in the 
Cartesian coordinate. In the other method, the invariant mass of the pair is constructed 
in the r-z plane with an opening angle cut of (0e+e- < tt/IO) in the azimuthal plane, here 
the z-axis is the beam direction. The random combinatorial background is constructed 
by rotating the partner candidate lepton by 180°. Raw yields of the photonic electrons 
are the excess over the combinatorial background in the low mass region (Me+e- < 150 
MeV). Comparing these two methods, the second one with an opening angle cut has 
the ability to reduce the combinatorial background and the effect of the secondary vertex 
resolution, which shift the mass to higher region (the second bump observed in the first 
method), see Fig. 13. 9[ 

From the measurement of the photonic invariant mass with the combinatorial back- 
ground subtracted in each px bin, the photonic electron raw yield is obtained. Electrons 
could be randomly combined with hadrons due to hadron contamination from the dE/dx 
selections. The different dE/dx cuts were tried to estimate the systematic uncertainties. 
The uncertainties of the combinatorial background (5 — 15%) dominate the systematic 
error for photonic electron spectra. The photonic background raw yields need correction 
by the efficiencies and detector acceptance as what we did for inclusive electron analysis. 
In addition, the partner finding efficiency, which was calculated by taking the ratio of 
e~^{e~) with partner e~{e~^) found to total converted electrons, was done from the same 
reconstruction process in embedding. Since the embedded 7r° was flat in pt, the charged 
71 Pt distribution measured by STAR |Lqe06| was used to weight the efficiency, shown in 
Fig. [3101 

From the same embedding data, EMC electron analysis gives the consistent back- 
ground efficiency, shown in the right panel of Fig. 13.111 At TOF electron p^ region, the 
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Figure 3.9: Left Column: Photonic invariant mass from TOF at low px- Right 
Column: Photonic invariant mass from TPC at high pt- Upper Raw: Photonic invariant 
mass reconstructed in the Cartesian coordinate. Bottom Raw: Photonic invariant mass 
reconstructed in the r-z plane with an opening angle cut of (0 < vr/lO) in the azimuthal 

plane. 
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Figure 3.10: Panel (a): Converted electron px distribution (blue solid line) and 
electron pt distributions with partner track found (red dashed line). Panel (b): 
Weighted electron pt distributions. Panel (c): Background reconstruction efficiency. 
Panel (d): Weighted background reconstruction efficiency. 
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effect of the momentum resolution to the background efficiency is smalL A random sam- 
phng method was also used to cross check the normal weighted method. The normalized 
Pt spectra shape of 7r+ was used to sample the electron candidates randomly in each pt 
bin. Both the two methods give consistent efficiencies, shown in the left panel of Fig. 13.111 
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Figure 3.11: Left panel: Comparison of different ways to derive background efficiency 
at TOF electron pt region from embedding. Right panel: Comparison of the 
background efficiency from different electron analysis. 



To systematically estimate the uncertainties of different tt spectra for the background 
efficiency weig hting, the PHENIX 7r° spectrum |PhpiO| was also used to compare with 
STAR TT"'" spectrum, shown in the Panel (c) of Fig. I3.12[ The difference of the efficiencies 
weighted by STAR Tr"*" and PHENIX 7r° pt distributions is small. At pt around 1-5 GeV/c, 
the efficiencies are nearly the same. Below 1 GeV/c, the difference is up to 20%. At high 
Pt (> 5 GeV/c), it is -5%, see the Panel (a),(b) of Fig. [SH 

Simply, the techniques are just shown in detail for minbias 0-80% and central 12% 
collisions, since good statistics can be collected in these two data sets. The techniques for 
other sub-centralities from minbias events are the same. 

After corrected by detector acceptance and efficiencies, the inclusive and photonic 
electron spectra were shown in the Panel (a) of Fig. 13.131 The ratios of the inclusive spectra 
and the photonic electron spectra are consistent with unit at low px (< 0.9 GeV/c). The 
increasing ratios above unit at pt > 0.9 GeV/c indicate that the observable signal, shown 
in the Panel (b). Due to poor statistics of the invariant mass by selecting tagged electrons 
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Figure 3.12: Panel (a): Comparison of tlie background efficiency from different 
weighting. Panel (b): The ratio of the efficiencies. Panel (c): Comparison of PHENIX 

7r° spectrum and STAR 7r+ spectrum. 
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from TOF, we use TPC dE/dx instead. Panel (c) shows the comparison of the photonic 
background spectra from TOF and TPC. Within errors, they are consistent. 
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Figure 3.13: Panel (a): Centrality dependence of inclusive electron and photonic 
electron spectra. Panel (b): The ratio of inclusive electron over photonic electron 
spectra. Panel (c): Comparison of photonic spectrum from TOF and that from TPC in 

central 12% collisions. 

Table 13.31 summarizes the cuts for tagged electron and its partner electron for the 
conversion pair selection. 

The final spectra were also corrected by the efficiencies of the ncTg cuts. The efficiencies 
of dctte+e-, Me+e- dJid Opening angle cuts were included in the background reconstruction 
efficiencies from embedding. From PYTHIA study, only ~2% electrons from 7r° and r] 
Dalitz decay are lost by cutting on Me+e- < 150 MeV, shown in Fig. I3.14[ Including 
vector meson decays, a 2-3% contribution from other sources to the total background was 
taken into account [STAcOS] . 

Since the tighter cuts for vertex-Z and eta could reduce the material influence, which 
is sensitive to the ratio of signal to photonic background, another cuts of \Vz\ < 20 cm 
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Figure 3.14: Invariant mass distributions of vr" and f] Dalitz decays from PYTHIA. 
The lost fraction above the mass cut is ~1.4% for all electron px- For electron pt > 1 

GeV/c, the lost fraction is ~2%. 
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Figure 3.15: Comparison of the inclusive and photonic electron spectra and their 

ratios with different Vz and t] cuts. 
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Table 3.3: Electron cuts for photonic background analysis 



Tagged electron 


partner electron 


inclusive electron cuts required 


nFitPts > 15 




nPitPts/nMax > 0.52 


-1 < riCTe < 3 (TOF) 


— 1 < nae < 3 


0<nae<2 (TPC) 




dca^+e- < 1 


cm 


Me+e- < 150 MeV 


Me+e-(e) < 150 MeV 


Opening angle (0) < vr/lO 



and —0.7 < \ri\ < were applied to cross check the final inclusive and photonic electron 
spectra, shown in Fig. 13.151 The inclusive and photonic electron spectra are both little 
lower than those from former cuts, due to the reduced yields of converted electrons. The 
inclusive/photonic ratio is similar as central top 12%, whose Vz distributions is narrow 
from -20 cm to 20 cm. At higher px {> 2 GeV/c), the ratio from former cuts is a little 
lower than the central top 12% and the minbias result from new cuts, which is probably 
because the decrease of material by applying the new cuts. 

3.1.5 Centrality dependence of non-photonic electron 

After the photonic background subtraction from inclusive electrons, the non-photonic 
electrons, which are the signal from heavy fiavor decay, are extracted. Due to low pt 
huge photonic background and high px low statistics, the non-photonic electrons were 
measured covering p^ from 0.9 GeV/c to 4 GeV/c in minbias Au+Au collisions and pt 
from 0.9 GeV/c to 5 GeV/c in central 12% Au+Au collisions. 

Fig. 13.161 shows the centrality dependence of the non-photonic electron spectra. The 
bin-by-bin errors dominated by the uncertainties of inclusive electron yield extraction, 
photonic combinatorial background reconstruction and photonic background efficiency 
etc., were propagated to the final non-photoinc electron spectra. The overall systematic 
errors, like the tracking efficiency uncertainty etc., were counted in the total systematic 
errors, which are shown as the shaded boxes. 
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Figure 3.16: Centrality dependence of the non-photonic electron spectra. 
Table 13.41 lists the systematic error contributions. 



Table 3.4: Systematic errors contributions to the electron spectra. 



Inclusive electron yield extraction 


5 — 15% (point-to-point) 


Combinatorial background uncertainties 


minbias: 6 — 15% (point-to-point) 
Central 5 — 10% (point-to-point) 


ncTe cut: (0,3)(-l,3) - (-1,3)(-1,3), (0.2)(-l,3) 


3-5% 


Background reconstruction efficiency 


0.3 <Pt< 0.9 GeV/c up to 20% (weight) 
Pt > 0.9 GeV/c 3 - 5% (weight, EMC) 


rj Dalitz, Kaon decay, vector meson decay 


2-3% 


TPC tracking efficiency 


5% 


Central 10% 12% 


5% 
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3.2 Non-photonic electron azimuthal anisotropy distributions 

3.2.1 Heavy flavor anisotropic flow 

In non-central nucleus-nucleus collisions, the anisotropy in the initial coordinate space 
is transferred into final momentum anisotropy. The anisotropic fiow is strongly sensitive 
to the partonic rescattering. The elliptic fiow parameter V2 from Fourier expansion of 
the azimuthal distributions is considered as good probe to extrapolate early pressure and 
density |Sor99] . The mass dependence of f 2 for identified light hadrons can be well repro- 
duced by a hydrodynamical model in the low pt region (p-p < 2 GeV/c), which indicates 
that the collective motion evolves in early partonic phase of the collisions [HouOlj . The 
Pt dependence of V2 scaled by the number of constituent quarks is observed to be univer- 
sal. The scaling behavior can be illuminated by a quark coalescence model |MV03] . But 
for heavy fiavor quark, due to its extremely heavy mass, it can acquire fiow only when 
light quarks punch it very frequently in a very dense medium. The thermal equilibrium 
is approached through sufficient interactions in the bulk matter. So the measurement of 
f 2 for heavy quarks is vital to test the light fiavor thermalization and partonic density in 
the early stage of heavy ion collisions. 

The single electron from heavy fiavor hadrons semileptonic decay was used to infer the 
production of heavy quarks. The angular correlation between charmed meson and electron 
from its semileptonic decay is significantly strong. That indicates charmed meson f 2 can 
be indirectly measured through single electron V2 |GKR04l IDES04] . The non-zero single 
electron V2 was observed by PHENIX detector at RHIC |SS07j . STAR detector provides 
a large acceptance with \(j)\ < 2tt, \r]\ < 1.5. The measurement of single electron V2 from 
STAR will systematically enhance our understanding of the partonic thermalization in 
the high dense matter. But it is still a challenge due to the huge photon conversion from 
the virulent material in the STAR detector. 
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3.2.2 Inclusive and photonic electron elliptic flow 

In this analysis, we try to develop a method to measure single electron V2 from the 
data taken with the STAR experiment during the y/s^=200 GeV Au+Au run in 2004. A 
total of 12 and 19 million 0-80% minimum bias Au+Au events were used for the inclusive 
electron identification from the TOF and photonic electron reconstruction from the TPC, 
respectively. The vertex Z is required to be less than 20 cm in order to reduce the material 
influence as well as the pseudorapidity (—0.7 < rj < 0). 

Inclusive electrons are separated from hadrons with more than 90% purity after ap- 
plying a cut of their dE/dx (0 < na^ < 3) in the TPC and a cut of |l//3 — 1| < 0.03 in 
the TOF. 

As the photonic electron identification in the previous analysis for electron spectra, in 
the cylinder with beam axis and azimuthal plane, the invariant mass of the e^e~ pairs 
is reconstructed with an opening angle cut less than vr/lO. The photonic electrons are 
identified from the previous invariant mass subtracted by the combinatorial background 
in a very low mass region (Me+e- < 15 MeV). Other cuts such as Mg+e- < 10 MeV, 
Me+e- < 20 MeV are also tried to systematically study the photonic electron V2- 
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Figure 3.17: Left panel: The event plane azimuthal angle distributions from TPC and 
TOF events. Right panel: Event plane resolution from random sub-event method. 



51 



An event plane method was used for electron V2 |PV98] . The event plane was re- 
constructed from the final particles' azimuths excluding electron candidates {nae < or 
nae > 3) to remove the auto-correlations. The acceptance and efficiency of the detec- 
tors in azimuth was corrected by compensating the azimuth to a fiat distribution with 
weights from day-by-day run. Additional pt weights were also applied to improve the 
event plane resolution. The second order harmonic azimuth angle \l/2 for event plane can 
be calculated from the Q vector, as Eq. I3.2|3.31 



^2 = (^arctan /2, 



< ^2 < VI" 



Q i^Qxy Qy 



Wi ■ cos(20i), ^ Wi ■ sin(20i) | 

\ i i / 



(3.2) 



(3.3) 



Here, Wi is the weight for each track with both the track azimuthal compensation and pr 
weight contributions. 
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Figure 3.18: Inclusive electron distributions from TOF. 

The left panel of Fig. 13.171 shows the event plane azimuthal angle distributions from 
both TPC and TOF events. A constant function fit to the distributions gives reasonable 
X^/ndf, which indicates that the event plane was reconstructed successfully. 
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The event plane resolution was calculated using the random sub-event method |P V98j . 
Each event was divided into two sub-events "a" and "b" randomly. Two event plane 
azimuthal angles \l/2 and \E'2 were reconstructed correspondingly. Then the event plane 
resolution r = (cos[2(\l'2 — ^^rp)]) can be calculated from Eq.(14) and (11) of |PV98] : 

(cos[2(vl/2 - ^rp)]} = 0^X2eM-xlm X [loixlm + hixl/^)] (3.4) 
(cos[2(v[/^ - XI/,,)]) = ^(cos[2(vI/i-vl>^)]) (3.5) 
X2 = v^/a = V2V2N (3.6) 

"^rp is the event plane angle calculated from the event plane Q vector. Iq and Ji 
are the Bessel functions. The event plane resolution is around 72% from random sub- 
event method, shown in the right panel of Fig. 13.171 The event plane angle "^rp is used 
as reference to obtain the azimuthal distributions {dN/dA(j)) of inclusive and photonic 
electrons. Fig. 13.181 shows the dN/ dA(f) distributions for inclusive electron in two typical 
Pt bins. The f 2 can be extracted from fitting the dN/ dAcj) distributions in each px bin 
corrected by the event resolution. 

The electron pairs in the mass window Me+e- < 15 MeV are filled into dN/dA(j)aii 
distributions as well as the combinatorial background electrons dN/dA(l)com- Then the 
photonic electron dN/dAcpp^ distributions were derived as: 



dN/dA(j)ph = dN/dA(Paii - dN/dA(j),om (3.7) 

Fig. 13.191 shows the extraction of the dN/ dAcf) distributions for photonic electrons. 
Photonic electron f 2 was also obtained from fitting to the dN/ dA(j) distributions in each 
Pt bin corrected by the event resolution. 

The inclusive and photonic electron V2 are shown in Fig. I3.2U[ They are higher than the 
meson and baryon f 2 shown in lines |DES04j in the low pt due to the decay kinematics. 
The different invariant mass cuts are systematically tried for the photonic electron V2- 
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Figure 3.19: Left panel: Extraction of photonic electron dN/dA(f) distribution from 
TPC. Right panel: Photonic electron dN/dAcf) distribution (zoom in the photonic 
electron A0 distribution in the left panel). 
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Figure 3.20: Inclusive electron V2 from TOP and photoinic electron V2 from TPC 
comparing with meson V2 and baryon V2- Error bars are only statistical. 
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The statistical errors are limited by inclusive electron yield. 



3.2.3 non-photonic electron elliptic flow 



The non-photonic electron f 2 and its propagated error can be derived from the formu- 



lae: 



pho 



non 
'^2 



r — 1 



VCT^.inc 

r — 1 



where r is the yield ratio of inclusive over photonic electrons, which is from the mea- 
surement of their spectra. The statistic error of photonic electron is extremely small to 
be neglected comparing to that of inclusive electron. The non-photonic electron V2 is 
shown in the bottom panel of Fig. I3.21[ The wild error bars are due to the very small 
r(~ 1.3 — 1.5 at 2 — 3 GeV/c) and large errors of inclusive electron V2- The small r is 
due to the huge photon conversion from the material in the STAR detector. The poor 
statistics of inclusive electron is due to the small acceptance of current TOF tray. 
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Figure 3.21: Current status of TOF non-photoinic electron V2 from inclusive electron 
V2 with the photonic background subtracted. 



The coming update for STAR detector will give us a very good chance to measure non- 
photonic electron V2 precisely. The hydrodynamic model and quark coalescence model for 
heavy flavors will be tested from this future experiment at STAR. 



55 



CHAPTER 4 



Charm energy-loss, freeze-out, flow properties and 

cross-section 

In this chapter, the measurements of Kn at low (< 2 GeV/c) in 200 GeV 

Au+Au colhsions will be introduced. In addition, we use a newly proposed technique 
to identify muons from charm decays at low pt [LZZXOG] . Combined with non-photonic 
electron measurement, all three measurements together stringently constrain the total 
charm production cross-section at mid-rapidity covering ~ 90% of the kinematics. They 
allow us to extract the charmed hadron spectral shape and to study the charm energy-loss 
and radial flow properties. 

4.1 D'^ reconstruction in Au+Au collisions 

A total of 13.3 million 0-80% minbias triggered Au+Au events in year 2004 Run IV 
were used for the (-0°) direct reconstruction through hadronic channel {D^ K~7c^, 
fC+vr"). The decay branching ratio is 3.83%. The collision vertex Z was required 
from —30 cm to 30 cm. In what follows, we imply {D^ + D^)/2 when using the term 
unless otherwise specified. 

Without the inner tracker devices, the exact decay vertex cannot be reconstructed 
due to insufficient track projection resolution close to the collision vertex. The invariant 
mass spectrum of mesons was obtained by pairing each oppositely charged kaon and 
pion candidate in the same event. The kaon and pion daughter tracks were identified by 
the dE/dx measured in the TPC The cuts for selection of the daughter candidates 
are listed in Table 14. 1[ 
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Table 4.1: Cuts for D'^ daughter candidates. 



momentum p 


> 0.3 GeV/c 


Pt 


> 0.2 GeV/c 


nFitPts 


> 15 


Global dca 


< 1.5 cm 


pseudorapidity r] 


(-1,1) 


pair rapidity 


(-1,1) 


ncr^ 


(-3,3) 




p < 0.7 GeV/c (-2,2) 




p > 0.7 GeV/c (-1,1) 



Table 14.11 lists the cuts for the kaon and pion candidate tracks. 




Figure 4.1: Left panel: Kaon/pion invariant mass distributions from same event (red 

open circle) and invariant mass distributions from mixed event (blue curve). Right 
panel: Kaon/pion invariant mass distributions from same event after the combinatorial 

background subtracted. 

The random combination of kaon and pion pairs contribute to a huge combinatorial 
background. An event-mixing method |Hbz03l ISTAc05j was provided to represent the 
random combinatorial background. In this method, the kaon in one event and the pion 
in another event from the event buffer were selected to generate a reference event-mixing 
invariant mass distribution. The buffered event candidates were required to have similar 
event environment. The mixing-buffer is divided into 10 RefMult bins and 10 Vertex Z 
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bins. Each RefMult bin and Vertex Z bin has roughly the same number of events. Each 
event is mixed with two other events in the same RefMuh bin and Vertex Z bin. Due 
to the muhi-combination with different event buffers, the statistics of the combinatorial 
background can be increased ~ 4 times. Then the event-mixing spectrum was normalized 
with a factor obtained by comparing the entries in the two spectra with invariant mass 
> 2.5 GeV/c^. In order to increase statistics, and harD^ signals are added together. 
Left panel of Fig. 14.11 shows kaon/pion invariant mass distributions from same events 
and mixed events after normalization. Due to the huge magnitude of the combinatorial 
background, the signal is invisible. But if the invariant mass distributions from same 
events are subtracted by the event-mixing combinatorial background, the signal can be 
seen within a mass region of 1.68—2.05 GeV, shown in the right panel of Fig. 14. 1[ The 
small peak near mass ~ 0.96 GeV is known as K* . 
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Figure 4.2: Extraction of the signal raw yields by fitting with open- width and 
fix-width from kaon/pion invariant mass distributions. 

After the event-mixing background subtraction, the kaon/pion invariant mass distri- 
butions within the mass window (1.68—2.05 GeV) were shown as the blue solid circles in 
Fig. 14.21 The red open circle is the invariant mass subtracted by the residual background. 
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A gauss+linear or second order polynomial (Pol-2) function as 



W X An -(m-m)2 

/(m) = X + (4.1) 

V27r X cr 

was used to fit the data to extract the signal, where W = 0.005 GeV is the bin width. 
Aq is the signal raw yield. /' is the linear or the Pol-2 function, a is the width of the 
gaussian mass distribution. From embedding, the width of the mass is around 13 MeV, 
but in real data, the width turns out to be smaller. Since in suh-px bins, due to statistical 
limit, the mass distributions can not be fitted by open width. So the width was fixed as 
10 MeV, and the overall uncertainty of fix-width and open-width fits, which is around 
10%, was done at < 3 GeV/c for good statistics, shown as Fig. 14.21 

To estimate the uncertainties of the combinatorial background, we use linear and Pol-2 
function to describe the background by varying several fit mass regions. The average yield 
from linear and Pol-2 fit with best was extracted as the final raw yield. The dominant 
bin-by-bin systematic error was estimated as the largest deviation from all the fits with 
reasonable x^? since all the fit results are considered as a uniform distribution. 

Fig. 14.31 gives the fits for three sub-p^ bins using linear function (the first row) and 
Pol-2 function (the second row) with best x^- 

Table 14.21 lists the systematical uncertainties for spectrum in minimum bias 200 
GeV Au-|-Au collisions. 



Table 4.2: D'^ systematical uncertainties in minimum bias 200 GeV Au+Au collisions. 



Pt 


bin-by-bin 


op en /fix- 


Full Field / 


Tracking 


Track cuts 


Total 


(GeV/c) 


fit 


width 


Reversed FF 


eff. 


(dE/dx) 


sys. 


0.2-0.7 


42.9% 










49.5% 


0.7-1.2 


26.8% 


10% 


9% 


10% 


18% 


36.4% 


1.2-2.0 


27.0% 










36.6% 



The final spectrum in 0-80% minbias AuAu collisions at 200 GeV are shown in 
Fig. HE 
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Figure 4.3: Kaon/pion invariant mass distributions by fitting with linear and Pol-2 

function in eacli px bin. 



60 



4.2 Muon from charm decay at low transverse momentum 

In recent experiments in STAR, due to large random combinatorial background, it is 
difficult to reconstruction charmed hadrons directly. And single electron measurement is 
impossible due to overwhelming photon conversions in the detector material and 7r° Dalitz 
decay at low pt, where the yield accounts for a large fraction of the total cross-section. 
Nevertheless, the charm total cross-sections have been measured in d-l-Au collisions at 
RHIC by a combination of the directly reconstructed low px Kit and the non- 

photonic electron spectra |STAc05j . and by electron spectra alone |PHECA^ IPHECPPj . 
Although the systematic and statistical errors are large, the result indicates a much larger 
charm yield than predicted by pQCD calculations [STAcOSl ICNVOSj . Since most of the 
measurements to date at RHIC are from indirect heavy-flavor semileptonic decays, it is 
therefore important to find novel approaches to improve the measurements and also study 
in detail how to extract the maximum information about the heavy-fiavor spectrum from 
its lepton spectrum. 

In this section, we propose a new method to extract the charm total cross-section 
by measuring muons from charmed hadron semileptonic decay at low px (e.g. 0.17 < 
Pt ^ 0.25 GeV/c). Since muons in this pt range are a very uniform sample of the whole 
charmed hadron spectrum, the inferred charm total cross-section is insensitive to the detail 
of the charm spectrum. Once the cross-section is determined, the electron spectrum at 
higher px can be used to sensitively infer the charmed hadron spectral shape. Left panel 
of Fig. 14.41 shows the charmed hadron px spectrum before and after requiring its decayed 
muons at 0.16 < pr < 0.26 GeV/c. The similarity of the spectral shape shows that the 
muon selection reasonably uniformly samples the entire charmed hadron spectrum. The 
muons in this Pt range sample 14% of the charmed hadron spectrum. The muon yield 
is about 1/70 of the charmed yield due to an additional 9.6% {c ^ I + anything) decay 
branching ratio. Right panel of Fig. 14.41 shows the dependence of the muon yield on (pt) 
for a fixed total charm yield. The yield is normalized to yields at n = 10 and (pr) = 1.3 
GeV/c, see the power- law function Eq. 14.51 We also note that the muon yield has a very 
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weak dependence on n which demonstrates that over a wide range in (pt) , the muon yield 
is within ±15%. This is in contrast to the large variation of the electron yield integrated 
above pt of 1.0 GeV/c, where a factor of 8 variation is seen in Fig. 14.41 
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Figure 4.4: Left panel: Charmed hadron {D) spectra {dN / dydpT) as a function of pt 

at midrapidity before (solid line) and after (dashed line) a muon selection of 
0.16 < Pt < 0.26 GeV/c and \yi\ < 0.5. The later was scaled up by a factor of 70. Right 
panel: Lepton yields relative to the fixed total charm cross- section as a function of {pt) 
for a charmed hadron transverse momentum spectrum. Solid line shows muon yields 
with a kinematics selection 0.16 < pt < 0.26 GeV/c and \yi\ < 0.5. Dashed line shows 

electron yields with pt > 1.0 GeV/c. 



A total of 7.8 million 0-80% minbias and 15 million 0-12% central triggered Au+Au 
events in year 2004 Run IV were used for the charm-decayed prompt muons analysis 
covering a pt range of 0.17-0.25 GeV/c at mid-rapidity (—1 < rj < 0). The collision 
vertex Z was required from —30 cm to 30 cm. The number of fit points was required 
above 24. 

The low Pt single muons were analyzed by combining the dE/dx measured in the 
TPC and = (p//5/7)^ from the TOF at STAR |Sha06j . A variable called no"^, which 
is calculated as Eq. 14. 4[ was used for the muon dE/dx selection. 
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A 



nan 



X In 



dE I dx^ 



TPC 



(4.4) 



^/NdE/dx dE / dx^B^^^^^i 

where A is the cahbration factor which only basically depends on the detector in- 
trinsic character, N^E/dx is the number of dE/dx fit points for the track reconstruction, 
dE/dx^PQ is the energy-loss of track x measured in TPC, dE / dx^^^j^^^i is the predicted 
value for track x from Bichsel calculation. 

The muons were selected as —3 < na^ < 0, — 3 < na^ < —0.5 for the two pt bins: 
0.17 <pt < 0.21 and 0.21 < pr < 0.25 GeV/c, respectively. 
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Figure 4.5: Panel (a): Particle mass squared distribution = {p/l3/'j)'^ from the 
TOF measurements. A clear muon mass peak is shown on the left and the primary pion 
candidates are shown at the right. The pion contamination was estimated by applying 
pion samples using dE/dx selection. Panel (b): Inclusive muon DCA distributions 
were obtained by subtracting the residual pion DCA in the muon mass window. Panel 
(c): The single muon raw yield was obtained from a fit to muon DCA distributions with 
the background DCA distributions combining the primary particle DCA distributions. 



After the muon dE/dx selection, a clean muon peak can be identified within a mass 
window of 0.008 < < 0.014 GeV^/c^ measured from the TOF, see Fig. (a). A di- 
gaussian function was used to fit the muon distributions (purple curve). The residual 
pions contributing to the muon mass window were estimated by a pion distributions 
with |n(Tvr| < 0.5 (green histogram). After dE/dx and selections, the muon yields 
in the mass window were counted in each bin of the DCA distributions, see Fig. 14.51 
(b). The residual pions (green dotted histogram) were subtracted statistically from the 
DCA distribution within the muon mass window applied to the dE/dx selected pion 
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sample [LZZXOG] . The systematic uncertainties (18 — 25%) are dominated by the pion 
contaminations, which were estimated using different dE/dx cuts for the residual pions. 
After the residual pion contamination subtracted, the inclusive muon DCA distributions 
were obtained (red open circle). 

The dominant background muons are from pion/kaon weak decays at low px- Other 
sources of background (p fi^fi^,!] •yfi'^fi', etc.) are found to be negligible from 
simulations, due to their very low yield at low pt- The measurement of rj/iT ratio shows 
around 0.5 at high pt (> 3 GeV/c), but drops very fast at low px, its contribution to 
the low Pt muons is very small. The background muon from pion/kaon weak decays 
were subtracted using the DCA distribution from HIJING simulation. The single muon 
raw yield was obtained from a fit to muon DCA distributions with the background DCA 
distributions combining the primary particle DCA distributions jLZZXOG] . see Fig. 14.51 
(c). The Pt distributions for muon invariant yields in 0-12% central and minbias Au+Au 
collisions are shown as open crosses and diamonds in Fig. 14. 8[ respectively. 

To understand the background shape, some relevant checks were performed. The muon 
DCA distribution from pion decay and that from kaon decay have similar shape, shown 
in the left panel of Fig. 14. 6[ the relative yield of pion and kaon in HIJING dose not affect 
the background DCA shape. Embedding results give consistent shape of the background 
DCA, see right panel of Fig. 14. 6[ 



4.3 Decay form factors 

In this section, we discuss the charmed hadron semileptonic decay form factor and 
its effect on the lepton spectrum. Fig. 14.71 shows the electron momentum spectra from 
charmed meson decays at rest generated using the Particle Data Table |PDG02] . PYTHIA |STAc05i 
IPHECAAl IPHECPPI [SjOT] , pQCD calculations |CNV05j and from the CLEO preliminary 



inclusive measurement |CLEO04] . The spectrum generated by the PDG is according to 
the form factor of charmed meson decays to pseudoscalar K+l + u, vector meson K* + l + i' 
and non-resonance (Kit) +1 + 1/ where the K* mass is used for the (Kit) system. The 
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Figure 4.6: Left panel: Background muon DCA distributions from pion decay (circle) 
and kaon decay (histogram). Right panel: Background muon DCA distributions from 
HIJING simulation (circle) and embedding (histogram). 

decay partial widths (F) of the three dominant decay channels are: 

1. K + I + u with pseudo scalar meson in final state (Z^^ B.R.=7.8%) 

dV 



oc 



pI 



dq^ (l-gVM*2)2 

where is the invariant mass of the virtual W ul, px is the momentum of the 
kaon, and M* = 0.189 GeV/c^ is a parameterization of the effective pole mass in 
the decay. 

2. [Kn) + 1 + 1/ with non-resonant Ktt in final state(il'^ 4.0%) 

We use the K* mass (0.892 GeV/c^) for the Kit invariant mass and the form factor 
is the same as in the decay to the pseudo scalar meson. 

3. K* + 1 + 1/ with vector meson in final state(Z)^ 5.5%) 



.^f n oc^[(l-cosgO^|/74g^)p + |(l + cosgO^|/74g^)P + |sin^g,|/7o( 
dq^d cos Oi 3 3 



where 9i is the decay angle between the lepton and the vector meson, pv is the 
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vector meson momentum, 



M + m 



and 



1 4 A/f^n^ 

Ho{q') = ^[(M^ - - g^)(M + m)A{q') - ip-^A,{q^)] 
2mq M + m 

where Ai^2,V take the form of 1/(1 - q^/Mfy) with = 2.5 GeV/c^, 
2.1 GeV/c2 and ry = V{<d)/Ai{Q) = 1.62 ± 0.08, r2 = A2(0)/Ai(0) = 0.83 ± 
0.05 |PDG02l[Pde5I] . 

For different charmed hadrons, we assume that the relative branching ratios among these 
three channels are the same, and their decay electron spectra are the same. The overall 
charmed hadron to electron branching ratio r(c e) /T{c anything) is 10.3% |PDG02j . 
There is a possible ~ 5% difference between electron and muon decays due to phase space 
which was not taken into account in this analysis. Electrons at high momentum are mainly 
from decay channel (1) K+l+u because the kaon is lighter than the K* and the form factor 
of the decay channel to a vector meson {K* + / + z/) favors a low momentum lepton and 
higher momentum neutrino. Since PYTHIA uses a simplified vector meson decay form 
factor [SjOl] , it tends to produce a softer electron spectrum. Both the parameterization 
by Cacciari |GNV05j and formulae from the PDG agree with CLEO's preliminary electron 
spectrum. In addition, we also find that although the charmed mesons {D^ and D^) from 
\E'(3770) decay have a momentum of 244 MeV/c only and without correction of final state 
radiation |CLEO04] . it affects slightly its subsequent electron spectrum. 

We use our electron momentum spectra and that of Cacciari to generate electron spec- 
tra from charmed decay at RHIC. A power-law function of the charmed hadron trans- 
verse momentum spectrum was obtained from minimum-bias Au+Au collisions GQM05| . 



Fig. 14. 71 shows the ratios of those electron spectra divided by the spectrum using PYTHIA 
decay form factors |Sj01 , BTAc05j . The slightly soft form factor of the charm semileptonic 
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Figure 4.7: Left panel: Electron momentum spectra from charmed meson decays at 

rest. Histogram is the one with form factor from PDG. The dashed line is that of 
Cacciari's parameterization. The dotted line is from the simplified vector-meson form 
factor in PYTHIA. The solid line is that from the PDG which takes a form factor from 
^^(3770)— * D e. Right panel: Charm-decay electron spectra for three different form 
factors divided by the spectrum using the PYTHIA form factor. The histogram 

represents the ratio using the PDG form factor. The solid line is Cacciari's 
parameterization, and the dashed line is from PDG and takes a form factor from 

^(3770)^ D e. See text for detail. 

decay in PYTHIA convoluted with a steeply falling charm spectrum produces an electron 
Pt spectrum in Au+Au collisions at RHIC, which can be lower than the correct one by 
up to a factor of 1.5 at high pjn. Part of the discrepancy between experimental results 
and PYTHIA in electron spectra [STAcOSi IPHECA"A] can be explained by the decay form 
factor. Taking and from \E'(3770) decay as if it were at rest, results in slight change 
on the electron spectrum. 

In the following, if we say 'form factor decay' that means we let the charmed hadrons 
decay using the form factor discussed above. 



4.4 Combined fit: Charm cross-section, freeze-out and fiow, 
energy- loss 

At RHIC energy, charm total cross-section have been measured from single electron 
spectra alone at pt > 0.8 GeV/c by PHENIX. But the charmed hadron pt spectrum is 



67 



unknown and the electron pt region is corresponding to higher px of the charmed hadron 
spectrum. Then the charmed hadron at low pt where the yield accounts for a large 
fraction of the total cross-section, is missing. The method to estimate the background from 
cocktail is also strongly model dependent. The direct extraction from electron spectrum 
will have large systematic uncertainties. STAR provided a method to extract the charm 
production cross-section by a combination of the directly reconstructed low px Kit 
and the non-photonic electron spectra |STAc05j . The measurement constrains most of 
the cross-section. But due to large combinatorial background, the systematic uncertainties 
are still large. 

On the other hand, as discussed in the previous sections, STAR collaboration has mea- 
sured D° hadronic decay channel (0.2 < pt < 2 GeV/c), single muon with high precision 
at low Pt (0.17 < pr < 0.25 GeV/c) and single electron (0.9 < pr < 5 GeV/c) at 200 
GeV Au+Au collisions. Charm cross-sections at mid-rapidity {da^^^ / dy) were extracted 
from a combination of the three measurements covering ~ 90% of the kinematics. In this 
section, the detail of how we perform the fit combining all the data points from these three 
measurements to extract the charm cross-section will be presented. In addition, the ex- 
traction of charmed hadron freeze-out temperature and flow velocity based on blast-wave 
model will be discussed. 

Firstly, the decay kinematics of charmed hadron to electrons were studied in pre- 
vious section. The same improved charm semileptonic decay form factors were used 
for all charmed hadrons. Assume similar pt spectrum shape between different charmed 
hadrons. And assume their decay electron spectra are similar. We applied the D° mass 
(1.863 GeV) for the form factor decay. The decay electron spectra were normalized by 
the fraction in total charmed hadrons from e~^e~ collisions at y/s= 91 GeV from 
PDG |PDG04 j: R = Noo/Ncc = 0.54 ± 0.05, and the charm branching ratio to electron 
(c ^ e B.R. = 10.3%). There is ~15% systematical difference between charm decay to 
electron and muon. In the following combined- fit, the single electrons are all from charm 
decays (Bottom contributions will be discussed in next chapter). 

A power-law function is used to create charmed hadron px spectra. The function takes 
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the form: 



dN rfiV2(n-l)(n-2),^ ^ 2pT 



27rdypTdpT dy n{n — ?>)'^{pt)'^ {PT){n — 3)' 
where dN/dy is the yield and n and (pr) are the parameters controUing the shape of 
the spectrum. A pt distribution with a set of these parameters was used as input 
for form factor decay, and the decay electron spectrum was obtained. A 3-dimensional 
scan on the {dN/dy, (pr), n) "plane" was done to fit and muon/electron data points 
simultaneously. The point with the smallest value was set to be the fit result. When 
we calculate the x^, we try to exclude the correlations between different measurements [?]. 
Correlated errors are, however, treated explicitly when there are a number of results of the 
form ± (Tj ± A that have identical systematic errors A, where cTj includes the statistical 
error and uncorrelated (bin-to-bin) systematical error. In this case, we use the quadratic 
sum of these errors (erf ± Af)-*^/^, where the modified systematic error 

A, = a,A[S(l/a|)]V2. (4.6) 
Then x^ was calculated from the following equation. 



X 



where yoi y^n ye denote the measured yields of D^, muons and electrons, cr^), cr^, Ue denote 
the measured errors, where we use statistical errors only to calculate the statistical errors 
for the fit results and use total errors calculated as above to estimate the systematical 
errors for the fit results, /d, fn, fe denote the expected values from input power law 
function for and its decay curve respectively. To avoid the pr position issue in large 
Pt bins, we used the integral yield dN instead of dN/pxdpx in each px bin. 

In addition, since the px distributions of D° is unknown, we also tried a blast- wave 
function as the input pr distribution. The blast- wave function is written as the 
following: 

dN mrcoshp pTsinhp 

oc / rdrmTKii — )Io[^f; ), (4.8) 



mTdrriT Jo Tfo Tfo 
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Table 4.3: The combined fit results for muons and electrons in Au+Au collisions. 



power- 
law 
fit 


measurements 




(Pt) (GeV) 


n 


X^/ndf 


mb (D^+fi+e) 


277 ± 26 ± 63 


0.92 ±0.06 ±0.12 


11.5 ±6.5 ±7 


18.6(4.2)/11 


0-12% (/i+e) 


311 ±26 ±64 


0.95 ±0.04 ±0.16 


13 ± 1.5 ±6.5 


18.5(0.6)/8 


blast- 
wave 
fit 


measurements 


da'yydy (fib) 


Tfo (MeV) 


(A) 


X^/ndf 


mb (D^+i^+e) 


271 ± 24 ± 57 


222 


0.27 (< 0.60) 


13.9(3.2)/8 


0-12% (/i+e) 


283 ±21 ±61 


220 


0.35 (< 0.63) 


45.2(2.3)/5 



where p = tanh^^jSt. The 3 free parameters are: dN/dy, freeze-out temperature Tfo and 
collective velocity Pm, where {Pt) = /5m x f • A 3-dimensional scan on the {dN/dy, Tjo, 
Prn) "plane" was also performed to fit D'^ and muon/electron data points simultaneously. 
The way to calculate errors and is the same as above. 

The error estimation was through the contour scan in the 3-D "plane" with the = 
Xmin + 1- The error of dN/dy was then obtained by projecting this 3-D contour into 
dN/dy axis. 

Once the dN/dy was extracted, the charm production cross-section per nucleon- 
nucleon interaction at mid-rapidity can be calculated from Eq. 14.91 



da, 



NN 
cc 



dy 



dN^io a, 



pp 



X Rx (4.9) 



In this equation, the factor R is the fraction in total charmed hadrons, as mentioned 
before. The number of binary collisions A^bin, which is from Glauber calculations, is 
293 ± 35 for 0-80% minbias Au±Au collisions and is 900 ± 71 for 0-12% central Au±Au 
collisions. 

Table 14.31 lists all the fitting results. 

Fig. 14. 81 shows , muons spectra and centrality dependence of non-photonic electrons 
spectra and the combining fit results for D^, muons and electrons spectra in Au±Au 
collisions. 

Within errors, both power-law fit and blast- wave fit give the similar da%^ /dy. There- 
fore, by averaged the two fits, da^^ /dy is presented to be 274±25(stat.)±60(sys.) /ib in 
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Figure 4.8: spectrum in minbias Au+Au collisions (star), muons spectra in minbias 
(diamond) and 0-12% (open cross) Au+Au collisions. Centrality dependence of 
non-photonic electron spectra. Solid curves show the power-law combined fit results. 
Dashed curves are from the blast-wave combined fit. 
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minbias Au+Au and 297±24±63 fxh in 0-12% central Au+Au collisions at ^/JNN=200 
GeV. The total charm cross-section per nucleon-nucleon collision (o"^^) following the 
method addressed in Ref. [STAcOSj is presented to be 1.40±0.11(stat.)±0.39(sys.) mb in 
0-12% central Au+Au and 1.29±0.12±0.36 mb in minbias Au+Au collisions at y'J/\r/v=200 
GeV. Fig. 14.91 shows the da^^^/dy as a function of A/'bin for minbias (i+Au, minbias 
Au+Au and 0-12% central Au+Au collisions. It can be observed that the charm cross- 
section seems to follow A^bin scaling from d+Au jSTAcOSj to Au+Au collisions which 
supports the conjecture that charm quarks are produced at early stages in relativistic 
heavy-ion collisions. However, the recent cross-section result from PHENIX is a factor of 
2 lower than STAR and the FONLL (NLO) calculations |CNV05l |Vog02i |Vog03| shown 
as the band, which under-predicts the minbias data by a factor of 5.1 ± 0.48(stat.) ± 
l.2{syst.yj;,\{theory). 
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Figure 4.9: Mid-rapidity charm cross-section per nucleon-nucleon collision as a 
function of number of binary collisions (A^bin) in (i+Au, minbias and 0-12% central 
Au+Au collisions. R factor is the fraction in total charmed hadrons. The solid line is 
from the average of the three values. Within the errors, the measured cross-sections are 
consistent with the number of binary collisions scaling. FONLL prediction, shown as the 
band, under-predicts the charm cross-section for collisions at RHIC. 
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The nuclear modification factor {RAuAu/dAu) |Adl02aj for single muon (open crosses) 
and non-photonic electron RAuAu/dAu (solid squares) are shown in Fig. l4.111 The RAuAu/dAu 
are obtained taking the ratio of the pt spectra in 0-12% central Au+Au collisions and the 
A^bin scaled decay electron spectra in d+Kn collisions. The decay d+Au electron spectra 
curve is from the combined fit with both TOF and EMC data [STACRAAl IPHECRAA] . 
see Fig. I4.10[ The bin-by-bin errors from the fit are propagated in the systematic errors 

of the RAuAu/dAu- 
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Figure 4.10: Single electron spectra measured from TOF are consistent with EMC 
results. In order to reduce the errors of the reference spectrum in d+Axi collisions, a 
combined fit was performed to fit both TOF data points and EMC data points. The 
variation of the power-law parameters in 1 — a was applied to propagate the 

uncertainties of the spectra shape. 

The muon RAuAu/dAu at low pt is consistent with unity considering uncertainties. The 
non-photonic electron RAuAu/dAu in 0-12% central Au+Au collisions is observed to be 
significantly below unity at 1 < < 5 GeV/c and is suppressed as strongly as that of 
light hadrons |Lqe06| , which indicates a large amount of energy- loss for heavy quarks in 
central Au+Au collisions. The measurement of non-photonic electron at high px from 
STAR EMC also shows strong suppression |STA(]RAA[ IPHEC^RAA] . Theoretical calcu- 
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lations [ArmOSt IDGW05] considering only the charm contributions to the non-photonic 
electrons agree with the measured non-photonic electron RauAu/clAu, while calculations 
with single electrons decayed from both bottom and charm quarks give larger RauAu/clAu 
values. The discrepancy from pQCD calculation assuming only gluon radiative energy- 
loss challenges our understanding of the detailed mechanisms of quark and gluon energy- 
loss in strongly interacting matter. Model calculations incorporating in-medium charm 
resonances/diffusion or coUisional dissociation can reasonably describe the non-photonic 
electron spectra [MT051 IHGR061 IRapOGj IVitOGj . see Fig. imi 
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Figure 4.11: Nuclear modification factor {RauAu/iIAu) of the spectra in — 12% 
Au+Au collisions divided by corresponding spectra in d+Au collisions. Bin to bin 
systematic errors are represented by the open boxes. Box on the right at unity shows 

the common normalization uncertainty in A^bin- Model calculations are presented: 
coalescence and fragmentation |Rap06] (double-solid lines), and coUisional dissociation 
of heavy meson |Vit06] (double-dashed lines). The solid, dashed and dot-dashed lines 
are Blast- Wave calculations with different freeze-out parameters of (T/o in MeV, {Pt)) 
for charmed hadrons (CH,1 — cr), multi-strange hadrons (MSH) and light hadrons (LH), 

respectively. 



The precise measurement of charmed hadron flow properties are expected to be a good 
test for partonic thermalization. The solid curve in Fig. 14.111 from bast- wave model with 
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the parameters (T/o=129 MeV, {Pt))= 0.477) for multi-strange hadrons cannot describe 
the data. The parameters obtained from pion, kaon and proton spectra (T/o=100 MeV, 
{Pt))= 0.6, dot-dashed curve) show even larger discrepancy [LZZXOGl IAda04a] . Due to 
the smearing of the charm semileptonic decay kinematics, only qualitative conclusion can 
be reached, that the charm spectra are not consistent with large flow and late freeze-out. 
Bast-wave parameters with low temperature and moderate radial fiow(dotted line), or 
with high temperature and low radial flow (dashed line) can describe our results. The 
data are also consistent with the dynamical models [MT051 IHGR061 Rap06 IVitOGj using 
finite charm interaction cross-section in a strongly interacting medium. This may connect 
the freeze-out parameters (temperature and flow velocity) to the drag constant in those 
dynamical models. Future upgrades with a direct reconstruction of charmed hadrons are 
crucial for more quantitative answers |STA05j . 

In summary, we report measurements of charmed hadron production at mid-rapidity 
from analysis of D ^ Kn , muons and electrons from charm semileptonic decays in 
minbias and central Au-|-Au collisions at RHIC. The transverse momentum spectra from 
non-photonic electrons are strongly suppressed at 1 < < 5 GeV /c in Au-|- Au collisions 
relative to those in d+An collisions. For electrons with px^2 GeV/c, corresponding 
to charmed hadrons with pj-r^ 3 — 5 GeV/c, the suppression is similar to that of light 
baryons and mesons. Detailed mo del- dependent analysis of the electron spectra with 
Pt < 2 GeV / c indicates that charmed hadrons have a different freeze-out pattern than the 
more copiously produced light hadrons. Charm differential cross-sections at mid-rapidity 
{da^^^ / dy) are extracted from a combination of the three measurements covering ~ 90% 
of the kinematics. The cross-sections are found to follow binary scaling as a signature of 
charm production exclusively at the initial impact. This supports the assumption that 
hard processes scale with binary interactions among initial nucleons and charm quarks 
can be used as a probe sensitive to the early dynamical stage of the system. 
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CHAPTER 5 



Discussion 

5.1 Heavy flavor energy loss 

As we discussed in the introduction chapter, the energy loss of heavy quarks is con- 
sidered as a unique tool to study the interactions between heavy quarks and the medium 
created in the heavy-ion collisions, and provide us important information of the medium 
properties. Experimentally, the high px modifications of heavy flavor hadron (D-mesons, 
B-mesons, etc.) yields are expected to be the direct variables to reveal heavy quark energy 
loss. But in current STAR experiment, it is very difficult to topologically reconstruct D- 
mesons or B-mesons. And due to large random combinatorial background, the same- and 
mixing-event method can only provide the px spectrum below 2 GeV/c with large sys- 
tematical uncertainties in minimum bias Au+Au collisions. In central Au+Au collisions, 
the combinatorial background becomes even larger, the signals are not good enough to 
extract useful information. Due to small acceptance, direct measurement of heavy flavor 
hadrons in PHENIX experiment becomes even harder. Nevertheless, the measurements 
of single electron from heavy quark decays were used as the indirect substitute. 

The single electron strong suppression, similar to light hadrons, was observed in recent 
measurements |STACRAA[ IPHECRAAj . Fig. 15.11 shows the nuclear modification factor 
of single muons (open crosses) and electrons (solid squares), RAuAu/dAu as a function 
of pt measured from STAR TPC+TOF detector. The RAuAu/dAu are obtained using 
d+Au fit curve as reference, which has been discussed in previous section. The PHENIX 
result IPHECRAAj and STAR EMC result [STACRAAj . Raa which use electron spectra 
measured in p+p collisions as references, are shown as open circles and solid circles for 
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comparison. Within errors, these results are consistent in different experiments. The 
RAuAu/dAu may show some Cronin effect |Cro73i IRua04bi IXinOS] due to enhancement in 
d+Au colhsions relative to p+p collisions RdAu > 1- But this effect is not significant 
compared to huge error bars. 
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Figure 5.1: Nuclear modification factors of single electrons as a function of pt- 



There are still some outstanding issues: how can the behavior of heavy flavor decay 
electrons reflect heavy quarks? Is the medium also opaque to heavy quarks? What is the 
contribution of bottom in single electron measurements? 

Fig. 15.21 shows the 2D-scattering plot of D-meson Pt(D) versus decay electron pr(e) 
from form factor decay. The inserted small panel shows the D-meson Pt(D) distributions 
when select electron at 1.45 < Pt{g) < 1-55 GeV/c. ~ 92% yields of D-meson are from 
1 < Pt{D) < 4 GeV/c. This indicates that the pt correlation between single electrons 
and their parent D-mesons is weak. 

On the other hand, due to large mass and small radiative angle, heavy quarks are 
predicted to lose less energy than light quarks via only gluon radiation - the " dead cone" 
effect |DK01j . Many theoretical calculations tried to explain the single electron strong 
suppression observed in experiments. A collisional (elastic) energy loss has been proposed 
to be taken into account for heavy quark energy loss calculations. 

Fig. 15.31 shows some recent theoretical calculations for single electron Raa from heavy 
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Figure 5.2: 2D-scattering plot of D-meson Pt(D) versus electron Pt{q) from form 

factor decay. 
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Figure 5.3: Theoretical calculations for single electron Raa compared to data. 
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flavor semi-leptonic decays. The DGLV radiative energy loss via few hard scatterings 
with initial gluon density dNg/dy = 1000 |DGV06j predicts significantly less suppression 
than data, shown as curve 1. The curve drops dramatically after including collisional 
energy loss |Wic05j . shown as curve 2, but still predicts less suppression than observed. 
Curve 3, which has a good agreement with data, is from the calculation of heavy quark 
energy loss via elastic scattering mediated by resonance excitations of D- and B- mesons 
and LO pQCD gluon exchange, and including heavy-light quark-coalescence at hadroniza- 
tion |Rap06] . The double curves indicate the uncertainties from the calculation assuming 
the D- and B- resonance width F = 0.4 — 0.75 GeV. Curve 4 shows the calculation of 
single electron suppression from collisional dissociation of heavy mesons in QGP by deriv- 
ing heavy meson survival and dissociation probability from the collisional broadening of 
their light cone wave function |Vit06j . The double-curve indicate the uncertainties from 
tuning the typical value of ~ 2 — 3. This calculation also shows consistent suppression 
with data. Curve 5 is for single electron only from D-meson decays with only radiative 
energy loss from BDMPS calculations via multiple soft collisions |Arm06] . Curve 6 is the 
same calculation as curve 5 except that it includes electron from B-meson decays. And 
assumes the transport coefficient g = 14 GeV^/fm |Arm06] . Both BDMPS calculations 
agree with data. 

There are still a few assumptions and uncertainties in theory. Several different pro- 
cesses can describe data within the large experimental uncertainties. The exact mechanism 
of heavy quark energy loss is still under intense theoretical and experimental investiga- 
tions. 



5.2 Heavy flavor collectivity 

Theoretical calculations have shown that interactions between the surrounding partons 
in the medium and heavy quarks could change the measurable kinematics [MTOSj IHGR06[ 
IDDZOGj . and could boost the radial and elliptic flow resulting in a different heavy quark 
Pt spectrum shape. Fig. 15.41 shows the picture that heavy quarks participate in collective 
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motion. We can simply treat heavy quark as an intruder, put into the hot medium with 
relatively very high density of light quarks. Due to their large mass and partonic density 
gradient, such a heavy quark may acquire flow from the sufficient interactions with the 
constituents of a dense medium in analog to Brownian motion. We expect that bottom 
quark has small collective velocity due to its extremely large mass and small interaction 
cross section. 




Figure 5.4: Heavy quarks (c, b) participate in the collective motion. 



Base on blast wave model |SSH93] discussed in previous section, a combined fit was 
performed to describe the spectral shape and to provide the freeze-out temperature T/o 
and collective velocity Fig. l5.5l shows the 2D distribution of Tfo versus Pm {{P) = f/^m) 
extracted from the combined fit to muon and electron spectra measured in — 12% central 
Au+Au collisions. The best fit with minimum gives the parameters as: Tfo = 220 MeV 
and (Pt) = 0.23. The 1-a (dotted curve) and 2-a (white solid curve) contours are from 
combined fits with statistical errors only. Black solid curve is for the l-cx contour from 
combined fit with both statistical errors and systematical errors. The fit results are not 
sensitive to Tfo, and {(3) is smaller than that of light hadrons. One can also see that Tfo 
and j3m are strongly correlated. 

To understand whether charmed hadrons may have similar radial flow as lighter 
hadrons, the blast wave fit result of charmed hadrons was compared to multi-strange 
hadrons and light hadrons (7r,K,p). Fig. 15.61 (a) shows these comparisons. Solid curve 
shows the blast wave fit for charmed hadrons (Tfo = 220 MeV and (Pt) = 0.23). The 
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Figure 5.5: Blast wave parameters Tfo versus (3m and 1-a contours from the combined 
fit to muon and electron spectra in 0-12% central Au+Au collisions. 

blast wave fit parameters for multi-strange hadrons (short-dashed curve) are Tf^ = 129 
MeV and = 0.48, and for light hadrons (short-dot-dashed curve) are Tfo = 100 MeV 
and (Pt) = 0.6. The long-dot-dashed curve, as a reference, is from a power law combined 
fit {{pj,) = 1.18 ± 0.02, n = 11.5 ± 0.5) to both TOF and EMC data in d+Au collisions, 
also shown in Fig. I4.10[ All these curves are scaled to match the measured cross sections. 
One may consider these comparisons highly model dependent. In principle, we can di- 
rectly compare the nuclear modification function between Omega and charmed hadrons 
since their masses are very similar. In practice, we take the blast-wave fit to Omega, 
which describes the Omega data well, and apply the charmed hadron semileptonic de- 
cay to obtain lepton spectra if the Omega and charmed hadrons would have the same 
freeze-out properties. This comparison is practically independent of model. 

The short-dashed curve in Fig. 15.61 (b) from blast-wave model with the parameters 
(Tfo = 129 MeV, {(3t) = 0.48) for multi-strange hadrons cannot describe the data |HYP06j . 
The parameters obtained from pion, kaon and proton spectra (Tfo = 100 MeV, (Pt) = 0.6, 
short-dot-dashed curve) show even larger discrepancy |Adl02aj . Due to the smearing of 
the charm semileptonic decay kinematics, only the qualitative conclusion that the charm 
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Figure 5.6: Panel (a): Different sets of blast wave parameters compared to measured 
lepton pt spectra. The solid, short-dashed and short-dot-dashed curves are blast-wave 

calculations with different freeze-out parameters of (Tjo in MeV, (A)) for charmed 
hadrons (CH), multi-strange hadrons (MSH) and light hadrons (LH), respectively. The 
long-dot-dashed curve is the power law decay curve fit to d+Kn data. Bin to bin 
systematic errors are represented by the gray bands for muons and open boxes for 
electrons. Panel (b): Nuclear modification factor {RAuAu/dAu) of the spectra in — 12% 
Au-|-Au collisions divided by corresponding power law decay curve fit to d+Kn data. 
Box on the right at unity shows the common normalization uncertainty in A^bin- The 
curves show the ratio of blast wave fit results and d+An power law, corresponding to 

panel (a). 
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spectra are not consistent with large flow and late freeze-out can be drawn. Bast-wave 
parameters with low temperature and moderate radial flow(dotted curve), or with high 
temperature and low radial flow (solid curve) can describe our results. This may indicate 
that charmed hadrons interact and decouple from the system differently from lighter 
hadrons. Future upgrades with a direct reconstruction of charmed hadrons are crucial for 
more quantitative answers |STA05j . 

5.3 Charm production cross-section — consistency and discrep- 
ancy 

We reported the charm cross sections from a combination of three independent mea- 
surements: D —>■ Kn , muons and electrons from charm semileptonic decays in min- 
bias and central Au+Au collisions at RHIC. While the discrepancy between STAR and 
PHENIX also exists. We perform a detailed comparison in this section. 

Charm total cross sections have been extracted from single electron spectra alone at 
Pt > 0.8 GeV/c by PHENIX. The single electron spectrum from charmed hadron decay 
peaks at pt ~ 0.5 GeV/c and therefore a large fraction of the total cross section is missing 
from the measurement. Recently, PHENIX reported new charm total cross section in p+p 
collisions from single electron spectrum at pt > 0.3 GeV/c. In all the cases, the inclusive 
electron spectrum contains photonic background. The photonic background level in this 
Pt range is comparable to STAR electron photonic background level at high pt- Cocktail 
and converter methods have been used by PHENIX to obtain the background. 

In this thesis, the charm cross sections are reported from combined three independent 
measurements: Kit, low Pt muons and single electrons in STAR experiments. 

Currently, the consistencies within STAR experiments are observed, while the discrepancy 
between STAR and PHENIX still exists. 

The comparison of the single muon/electron spectra measured from STAR TOP de- 
tector and the single electron spectrum measured from EMC detector in 0-12% and 0-5% 
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central Au+Au collisions, respectively, is shown in Fig. 14.101 These data points are scaled 
by A^bin since they are from different centralities. In order to see the comparison more 
clearly, we take the ratio of the data and the fit curve, see Fig. 15. 7[ Since EMC results 
have smaller errors, so it constrains the fit curve, thus the EMC data points are around 
unity and TOF result is a little higher, probably due to 0-5% has more suppression than 
0-12%. But within errors, the results from two independent measurements are consistent. 
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Figure 5.7: The comparison of the single muon (open crosses) and electron (open 
squares) spectra measured from STAR TOF detector and the single electron spectrum 
(open circles) measured from EMC detector in 0-12% and 0-5% central Au+Au 
collisions, respectively. The curve shows the fit to all the data points. Bottom panel 

shows the ratio of data and the fit curve. 

In Fig. 15.81 (a), the single electron spectrum from PHENIX measurement in p+p colli- 
sions, shown as triangles, has similar shape but systematically lower than STAR measure- 
ments, while STAR TOF result (circles) is consistent with STAR EMC result (diamonds) 
within errors. STAR TOF muon measurement (crosses) has also good agreement with 
STAR TOF electron measurement (squares) in 0-12% central Au+Au collisions. A form 
factor decay combined fit (solid curve) describes both low pt muons and higher pt elec- 
trons and gives the charm cross section at mid-rapidity as 297 mb. In addition, STAR 
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muons can not be fitted together with PHENIX electrons (scaled by A^bin) , see the dashed 
curve, which is mostly constrained by PHENIX electron data points with small errors. 
PHENIX published charm cross section is 123 mb. The cross section from fit to PHENIX 
data only is 131 mb. And the fit combining STAR Au+Au muons and PHENIX electrons 
gives 137 mb with large x^- The PHENIX cross sections are systematically a factor of 2 
lower than STAR. But both STAR and PHENIX are self-consistent in different collisions. 
Both observed the number of binary collisions scaling behavior of charm cross section. 
Fig. 15.81 (b), which shows the charm cross section at mid-rapidity as a function of number 
of binary collisions for STAR (circles), PHENIX (triangles) and FONULL calculations 
(solid line and dashed box), clearly summarizes the consistency and discrepancy among 
different data sets and theories. 
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Figure 5.8: Panel (a) Spectra comparisons between STAR and PHENIX. Spectra 
measured in p+p collisions are scaled by A^bin to compare with central results. Curves, 
from power law charmed hadron distributions via form factor decay, are used to fit data. 

Panel (b): The comparisons of charm cross section at mid-rapidity in different 
centralities and different measurements. FONULL calculation is shown as the solid line 
for comparison, and the dashed box stands for the uncertainties of theory. 



The FONLL and NLO pQCD calculations |(:;NVn51 |Vog02i |Vog03| underpredict the 
STAR minbias data by a factor of 5.1 ± 0A8{stat.) ± 1.2{syst.)^l'^{theory). As discussed 
previously [Xin05j . theoretical calculations tuned parameters to match the low energy 
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measurements and then extrapolated to high energies, and the parameters are probably 
energy dependent but not understood well yet. The charm cross section will be precisely 
measured with a direct reconstruction of charmed hadrons in the future, more detail will 
be discussed in the outlook chapter. 

The mid-rapidity charm cross-section per nucleon-nucleon collision {da^^^ / dy) can be 
converted to the total charm cross-section per nucleon-nucleon collision (o"^^) following 
the method addressed in Ref. [?]. The cr^^ is presented to be 1.40±0.11(stai(:.)±0.39(s?/s.) 
mb in — 12% central Au+Au and 1.29 ± 0.12 ± 0.36 mb in minbias Au+Au collisions at 
■s/sJfN = 200 GeV. But the different widths of the rapidity distributions from theoretical 
models will lead to substantial systematical uncertainties. Fig. 15.91 shows the charm cross 
sections as a function of rapidity distributions measured from STAR and PHENIX, com- 
pared to theoretical models |Vog04[ [RP031 [BCHOSj . Here the systematical uncertainties 
are dominant. PHENIX forward muon measurement [PHEMu] with large errors, shown 
as the solid triangle at rapidity ~ 1.6, gives consistent result with STAR, which is also 
significant higher than theoretical predictions. 
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Figure 5.9: Charm cross sections as a function of rapidity distributions measured from 
STAR and PHENIX, compared to theoretical models. 
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5.4 Identify bottom contribution in non-photonic electron spec- 
tra and V2 from Au+Au collisions at RHIC 



In the above discussions, the bottom contribution to the non-photonic electron spec- 
trum is neglected. We have discussed in the introduction chapter that the separation 
of bottom and charm contributions in current non-photonic electron measurements is 
very difficult. There are large uncertainties in the model predictions for charm and bot- 
tom production in high-energy nuclear collisions. And many theoretical calculations with 
charm or with both charm and bottom contributions describe the non-photonic electron 
suppression |DGV06l IWicOSl IRapOGj IVitOGl lArmOGj and V2 |GKR04l IZCK051 IHGROGj . 
Thus identification of bottom from the non-photonic electron measurements is crucial to 
better understand charm physics. In this section, we try a fit to non-photonic electron 
spectrum and estimate the bottom contributions. We also compare the V2 distribution 
from simulation to the experimental data and estimate the possible charm V2- 



5.4.1 Fit to non-photonic electron spectrum and relative cross section ratio 

The non-photonic electron spectrum up to 10 GeV/c has been measured by STAR 
experiment in 200 GeV p+p collisions. The idea is that we use the sum of electron spec- 
tra from both charm and bottom decays in PYTHIA model |Sj01] to fit the STAR p+p 
data [STAGRAA] to extract the fraction of the bottom contribution. As we discussed 
previously, the D-mesons and their decay electrons spectra from default PYTHIA param- 
eters are soft [LZZX06] . see Fig. 14. 7[ Thus a modified Peterson Fragment Function (FF) 
and the high pt tuned parameter are used to make spectra harder to be comparable with 
the form factor decays |LXY04] . 

Table EH] lists the parameter initialization for PYTHIA 6.131: 

Fig. 15.101 (a) shows the px distributions of the heavy fiavor hadrons and their decay 
electrons from PYTHIA with above parameters. The D-meson spectrum, shown as the 
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Table 5.1: PYTHIA parameters for heavy flavor decays. 



Parameter 


Value 


iViOil/ij 


4 (charm), 5 (bottom) 


quark mass 


= 1.25, mt = 4.8 (GeV) 


parton dist. function 


CTEQ5L 


scale 


4 


K factor 


3.5 




1.5 


Peterson Frag, function 


e = 10-^ 


high -pT tuned PARP(67) 


4 



hatched band, is normalized to 

dN/dy = dN/dy{D^)l{Nun)lRdAulR, (5.1) 

where dN/dy{D°) = 0.028±0.004±0.008 measured in d+Au coUisions |STAcn5] . (Nun) = 
7.5 ± 0.4 in d+Au collisions. RdAu = 1-3 ± 0.3 |Xin05] . R factor stands for fraction 
in total charmed hadrons, the fragmentation ratio R{c —>■ D^) = N^o/Ncc = 0.54 ± 
0.05 |PDG04j . All these normalization errors are propagated into the uncertainty band of 
the D-meson spectrum. The curve in this band is the lower limit of the D-meson spectrum 
in our simulation. Correspondingly, its decay electron spectrum is shown as the solid band. 
The non-photonic electron spectrum measured in p+p collisions at STAR [STACRAA] 
is shown as the open squares. The decay electron band alone can describe the data, 
indicating that the contribution of electrons from bottom decay could be very small. In 
order to estimate the upper limit of bottom contribution, we use the lower limit of the 
decay electron spectrum, shown as the open circles. B-meson spectrum (solid curve) and 
its decay electron spectrum (open triangles) are normalized by varying the ratio of (j),i/acc- 
The summed spectrum (solid circles) by combining the lower limit of D ^ e and i? ^ e is 
used to fit STAR data in p+p collisions, and then the upper limit of S — e contribution 
will be extracted. 

Fig. 15.111 (a) shows the fit ^ function of the unique variable a^jocc- The best 

fit with a minimum x^/ndf = 16.6/14 gives the upper limit of the total cross section 
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Figure 5.10: Panel (a): D-/B- mesons and their decay electron spectra from PYTHIA. 
The B + D ^ e fit to STAR non-photonic electron data in p+p collisions. Panel (b): 
The relative spectra ratio, upper limit of -B — > e contributions as a function of pt- 

ratio as af,i/acc < (0.49 ± 0.09 ± 0.09)%. The first term of the errors is calculated from 
= Xmin + 1- The second term is from the 15% normalization error of the dN/dy 
converted to total cross sections due to the uncertainties of the model dependent rapidity 
distributions [Xin05j . Fig. 15.111 (b) shows the B- /D- mesons rapidity distributions from 
PYTHIA. The cross section ratio from FONLL calculation is 0.18%-2.6% |CNV05j . The 
upper limit is consistent with theory prediction. 

The upper limit of i? — e contributions as a function of pt is shown in Fig. 15.101 (b). It 
is increasing and becomes flat around 7 GeV/c. The pt crossing point, where the bottom 
contribution is equal to charm, of electron spectra from B,D decay is very sensitive to the 
cross section ratio, since at high px-, these electron spectra shapes are similar. From the 
B + D ^ e fit to STAR p+p data, we estimate the crossing point p^>7 GeV/c. 

Table 15.21 lists the crossing points of heavy flavor decay electrons in several Pt bins. 
Table 5.2: Crossing points of heavy flavor decay electrons as a function of Pt- 



Pt (GeV/c) 


2 


3 


4 


5 


6 




8 


{B^e)/{D^e)< 


0.11 


0.31 


0.53 


0.77 


0.85 


1.2 


1.1 
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Figure 5.11: Panel (a): Fit a function of (Thi/acc- Straight lines is for the 

= Xmin + 1- Panel (b): B- (solid curve) /D- (dashed curve) mesons rapidity 

distributions from PYTHIA. 

5.4.2 Fit to non-photonic electron V2 

Besides the non-photonic electron spectrum, the non-photonic electron V2 has also 
been measured in 200 GeV Au+Au collisions at RHIC [?]. In this measurement, bottom 
contribution has not been separated, which can be studied by comparing simulations and 
data. Since heavy flavor hadrons pt distributions and V2 are unknown, our simulations 
have to base on the following assumptions: 

— The same relative (-B — e)/(D — e) ratio from p+p to Au+Au. 

— Assume the B-/D- meson V2 as the inputs for the simulation, here we assume three 

aspects: 

• I: B-/D- meson V2 are similar as light meson V2- 

• II: D-meson V2 as light meson V2 but B-meson does not flow. 

• III: B ^ e contribution is neglected and D-meson V2 decreases at pt > 2 
GeV/c. 
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Here heavy flavor baryons, Ac, are taken into account as 10% of total lieavy fla- 
vor liadrons |PDG04[ ?]. Tlieir V2 are assumed to follow liglit baryon V2. Tliis baryon 
contribution effect in tfiis simulation is small. 



We use the light meson V2 curve from fltting experimental data [DES04j as the input 
B/D V2 distributions (Assumption I), see Fig. 15.121 (a). That means in each pt bin, the 
B/D A</) distribution is initialized. The electron A0 distributions in each px bin will be 
obtained via B/D decays in PYTHIA model. Then the electron f2, shown in Fig. 15.121 
(b), will be extracted by fltting the distributions in each pt bin. 
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Figure 5.12: Panel (a): Assumed B-meson V2 (open circles) and D-meson V2 (dashed 
curve) as light meson v'^. Panel (b): Electron V2 from B-meson decays (open circles) and 

D-meson decays (open squares). 



Fig. 15.121 shows the obvious mass effect: The B/D V2 are assumed as the same, but the 
decay electron V2 can be very different due to decay kinematics |YFHP06j . This is not 
surprising, since we know B-meson is much heavier than D-meson and light hadrons. The 
decay electrons can only have a small momentum fraction of B-mesons. The momentum 
and angular correlations between decay electrons and B-mesons are weak, especially at 
low Pt- Therefore, at low pt the decay electron cj) angle will almost randomly distribute. 
So we see the zero or negative f 2 for the electron from B-meson decays. But from previous 
study, we know that bottom contribution below 3 GeV/c is small, thus the mass effect to 
the total electron V2 is not signiflcant. 
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Fig. 15.131 (a) shows the total electron f 2 from PYTHIA simulation compared to data. 
The measured non-photonic electron V2 from PHENIX is shown as the triangles. The 
solid curve (Assumption I) is the sum V2 of the two decay electron V2 distributions in 
Fig. 15.121 (b) by taking the relative ratio of {B —>■ e) / {D e) into account. It can not 
describe the data. If we assume B- meson does not flow (Assumption II), the total decay 
electron f 2 will become decreasing, shown as the band. The band is corresponding to the 
<^bb/^cc = (0.3 — 0.7)% (The upper limit, 0.49%, is in between). It has better agreement 
with data, but still higher. The decreasing of non-photonic electron f 2 could be due to 
5 — > e contribution and B- meson V2 could be very small. But below 3 GeV/c, i? — > e 
contribution is not significant. That indicates D- meson V2 should be smaller than light 
meson t>2 and start decreasing at higher p-p (> 2 GeV/c). 
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Figure 5.13: Panel (a): The total electron V2 from PYTHIA simulation assuming that 
bottom flows (solid curve) and bottom does not flow (band) compared to data. Panel 
(b): The total electron V2 from PYTHIA simulation fit to data and the estimated 

D-meson V2- 



So ignoring B ^ e contribution, we try to speculate the D-meson f 2 by fitting the data 
using decay electron t>2 (Assumption HI). In Fig. 15.131 (b). the best fit of the decay electron 
V2 is shown as the open circles. The estimated D-meson V2 is shown as the dashed curve, 
which is smaller than light meson V2 above 1 GeV/c and start decreasing above 2 GeV/c. 
This most possible D-meson V2 distribution shows that at p-r < 3 GeV/c, where the 
bottom contribution is negligible, D-meson has large ^2, indicating that charm strongly 
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flows in high dense medium, which could be the evidence of light flavor thermalization in 
QGP created at RHIC energy. 
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CHAPTER 6 



Detector upgrades and outlook 

6.1 STAR detector upgrades: TOF and HFT 

The experimental physics resuhs are inevitably limited by the detector resolutions, 
luminosities, acceptances and electronics etc. Physics developments require the upgrade 
of the detectors and electronics. The STAR Collaboration has proposed some important 
sub-detector upgrades. The full barrel Time-Of-Flight (TOF) detector [STA04] and the 
Heavy Flavor Tracker (HFT) [STA05J are two of them. The other detector upgrades like 
the Intermediate Silicon Tracker (1ST) and the Forward GEM Tracker (FGT) etc, will 
not be introduced in this thesis. 

The barrel TOF detector based on recently developed Multi-gap Resistive Plate Cham- 
ber (MRPC) technology will surround the outer edge of the TPC, covering — 1 < r/ < 1 
and ~ 2-7? in azimuth. The TOF system can achieve the required timing resolution < 100 
ps and the particle detecting efficiency > 95%. This detector will significantly extend the 
reach of the STAR scientific program, doubling the percentage of kaons and protons for 
which particle identification is possible to more than 95% of all those produced within 
the MRPC-TOF acceptance. Combined with existing STAR detectors, the barrel TOF 
detector will allow STAR to extract the maximum amount of information and statistics 
for heavy fiavor, dilepton and resonances physics, such as the single electron, muon 
spectra, the single electron t>2, J/"^, p, uj, vector mesons, etc. 

When combined with a possible future vertex detector upgrade, the proposed TOF 
detector will also reduce the integrated luminosity needed to measure a statistically robust 
sample of D^, D^, mesons by approximately an order of magnitude, enabling STAR 
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to make systematic studies of charm thermalization and meson flow. 

In addition, with the barrel TOF detector, it becomes possible to separate the bottom 
contribution in non-photonic electron measurement experimentally and charm physics will 
be more clearly understood. 

The expectant vertex detector HFT brings extremely high precision tracking capabil- 
ities to STAR with a resolution of 10 fim at the first layer of the detector, over a large 
pseudo-rapidity range, and with complete azimuthal angular coverage $ = 27r. The two 
layers of the Heavy Flavor Tracker (HFT), which are composed of monolithic CMOS pixel 
detectors using 30 /im. x SO/im square pixels, will be placed closed to beam pipe in STAR 
at radii of 2.5 cm and 7 cm, respectively. 

The HFT is designed to measure displaced vertices that are displaced 100 micros, or 
less, from the primary vertex. Therefore, the neutral and charged particles with very short 
lifetimes can be distinguished from the primary particles which originate at the collision 
vertex. The addition of the HFT will extend STAR's unique capabilities even further by 
providing particle identification for hadrons containing charm and beauty and electrons 
decaying from charm and beauty hadrons. Thus, the HFT is the enabling technology for 
making direct charm and beauty measurements at STAR. 

6.2 Direct measurements of charmed hadrons with HFT 

The Monte Carlo events in Au-|-Au coUisions at -y/s^ = 200 GeV are generated by 
using HIJING model. The standard GSTAR simulation package with modifications of 
the new detector design is used for the material geometry configurations. The new STAR 
tracking software package ITTF with HFT layers added is used for tracking. 

The physical interactiona between particles and the material of the detector are sim- 
ulated using the STAR implementation of the GEANT simulation package. This package 
is used in STAR and is a standard analysis tool which includes a detailed understanding 
of the TPC response function; including dead areas and reahstic detector resolutions and 
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responses. Realistic detector resolutions are used to smear the perfect position informa- 
tion, and the resulting simulated hits are used in tracking. For the SSD, the hits were 
smeared by 20 fim in and 750 fim in z, in agreement with the SSD specifications. For the 
HFT, the hits were smeared by 6 fim in both $ and z. 

In this section, we will focus on the charmed hadron reconstruction. The optimization 
of the HFT hit finding accuracy, vertex reconstruction performance, single track efficiency 
and ghost rate etc, will not be discussed. The details can be found in Ref. |STA05] . 

To demonstrate the power of the HFT, we have simulated several specific charm hadron 
decay channels, including K^tt^ and K^-k^-k^. Table W7[\ displays some of 

the properties of these channels. 



Table 6.1: Open charm hadron properties. 



Particle 


Daughters 


Cr (fim) 


Mass (GeV) 




K~'K+ (3.8%) 


122.9 


1.8645 


D+ 


ii'-7r+7r+ (9.5%) 


311.8 


1.8693 


Dt 


$7r+ (4.4%) 

TT+TT+TT- (1.2%) 


149.9 


1.9682 


Ac 


pK-Ti+ (5.0%) 


59.9 


2.2865 



Signal and background events are generated separately. The signal consists of one 
or per event. The transverse momentum distribution of the charmed hadrons follows 
a Boltzman distribution which reproduces the (pt) of D-mesons as measured by STAR 
in d+Au collisions at ^/s^ = 200 GeV |STAc05] and the rapidity distribution suggested 
by perturbative QCD calculations applying the program code Pythia |Sj01| . The back- 
ground is simulated using the MevSim event generator parameterized to reproduce the 
experimentally measured particle multiplicities in Au+Au collisions at ^/s^ = 200 GeV. 
Our parameterization is accurate for particles below 3 GeV/c. It may underestimate 
the background above this momentum. The distributions of reconstructed D-meson sig- 
nal and background were scaled to match the expected D-meson production per central 
Au+Au collision. Also, the higher track reconstruction efficiency in single (signal) events 
compared to central Au+Au collisions (background) was taken into account. 
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The heavy flavor tracker is designed to allow us to directly reconstruct mesons con- 
taining charm quarks. If this is true, we should see differences between the charm meson 
daughter tracks and the background and the primary tracks in several important vari- 
ables. The DCA between the tracks and the vertex is an important example of this; 
see Fig. 16.11 The distributions of reconstructed and daughters is compared to the 
primary track background. The charm meson tracks clearly have a broader distribution, 
driven by the decay length of the charm mesons. Cutting on the track DCA, then, will 
improve the signal to noise in the analyses. 

^ — I — I — I — \ — ^ — \ — I — I — I — ^ — \ — \ — I — I — I — \ — \ — \ — I — I — I — \ — \ — p 
refit with lielix + HFT liits 




dca (|im) 

Figure 6.1: The DCA distributions for (triangles) and (circles) reconstructed 
by HFT. The results of the dca distribution for primary tracks are shown in squares for 

comparison. 

The following selection criteria were also used to separate the D-meson signal from 
background. 

• The decay length 1: the distance between the primary vertex and the D-meson decay 
vertex (VO). 

• The DCA between the daughter tracks and the reconstructed decay vertex. 
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• Momentum angle correlation cos(0): 9 is defined as the angle between the D-meson 
momentum (vector sum of the two daughter momenta) and the vector joining the 
primary vertex to the D-meson decay vertex. 

The topological reconstruction of i^'^vr^'^Tr^ is introduced here as an example. 

The better signal to background can be obtained for the reconstruction of D*^ K^tt'^, 
since it is two-body decay. The HFT is designed as close to collision vertex as possible, 
so that it gives good primary vertex performance (~ lO/zm), which is good enough for 
the secondary vertex reconstruction. And the reconstructed TPC helixes are required to 
hit on the two layers of HFT. After refit the helixes with two HFT hits, the new helixes 
were obtained. The new helixes were used to calculate those variables mentioned above. 
Since it is three-body decay, decay vertex is estimated by calculating the average of 
three displaced points of K^tt^ , K^tt2 and ttj^ttJ. 

The reconstruction of charmed baryon A^ is similar as D^. But due to the smaller c,- 
and different charm fragmentation ratios, the cuts for A^, and should be different. 
The cuts are optimized to obtain good signal to background ratio. Table [6^ lists the cuts 
for the charmed hadrons. 

Table 6.2: Topological cuts for open charm hadrons. 



Cuts 


DO 


D+ 


Ac 


nFitPts > 


15 


15 


15 


\v\ < 


1.0 


1.0 


1.0 


DC A (global) > 




100 fiui 


35 fim 


DCA (VO) < 


35 fim 


100 fim 


40 fim 


decay length > 


150 fim 


150 fim 


50 fim 


cos(6') > 


0.996 


0.85 


0.92 



The dca, cos(^), and decay-length distributions for both primary tracks (open-squares) 
and decayed tracks are shown in Fig. 16.21 Clearly the decayed tracks are well separated 
from the primary track. For reconstruction, a slightly different method was used 
compared to that of the D^. The 'signal event' and 'background event' were mixed 
together. The distribution of the invariant mass from a K and two vr tracks is then 
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refit with TPC+SSD+HFT 



o D* decay tracks 




dca (|am) cos(e) decay-length (|am) 



Figure 6.2: The daughter tracks decay-length, dca, and cos(6') distributions 
(circles). The same distributions for primary tracks are shown as squares. 

formed. The number of tracks used in the background events is consistent with the top 
10% central Au+Au collisions at RHIC. The resulting invariant mass distributions for 
several pt bins are shown in Fig. 16.31 For the pt bins studied so far, the significance 
S/ \^S + B is better than 3. 



6.3 Bottom separation in heavy flavor measurements 

With the barrel TOF detector and vertex detector HFT, the D-meson px spectrum and 
non-photonic electron px spectrum will be measured precisely. The D-meson semileptonic 
decayed electron spectrum will be known based on models via decay kinematics. Compar- 
ison between the D-meson decayed electron and non-photonic electron px distributions 
will tell us the information of the bottom contribution in non-photonic electron measure- 
ments. This will also be tested from the measurements of D-meson f 2 and non-photonic 
electron V2 when the TOF detector and HFT are available. 

But before the barrel TOF detector and HFT installation in STAR, the separation of 
B-meson and D-meson due to their different life time (decay length) is proposed to be a 
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Figure 6.3: Invariant mass distributions for for several pt windows. The lines are a 

polynomial (up to 2nd order) + Gaussian fit. 

good try. Fig. 16 .41 shows the different decay length distributions of D-meson and B-meson. 

D-meson decay vertex will be closer to the collision vertex than B-meson. Currently, 
the existing vertex detector SVT and intermediate detector SSD will be used to help the 
reconstruction of the B-meson decay vertex. 
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Decay Length (i^m) Decay Length (i^m) 



Figure 6.4: Panel (a): D-/B- mesons decay length distributions. Panel (b): The 
fraction of B-meson in the total decay length of (B+D)-meson. 
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USA, 03/25/2006 - 04/01/2006. 

• Charm Cross Sections and Collectivity from its Semileptonic Decay at RHIC 

2nd International Conference on Hard and Electromagnetic Probes of High-Energy 
Nuclear CoUisions, Asilomar, California, USA, 06/09/2006 - 06/16/2006. 

• Identify Bottom Contribution in Non-photonic Electron Spectra and V2 from p+p 
and Au+Au Collisions at RHIC 

23rd Winter Workshop on nuclear Dynamics, Big Sky, Montana, USA, 02/11/2007 
- 02/18/2007. 



102 



Publication List 



• Y. F. Zhang (for STAR collaboration), "Open charm production in ^/lsNN) — 200 
GeV Au+Au collisions.", J. Phys. G: Nucl. Part. Phys., 32, S529-S532, 2006. 

• Z. Xu et al, "Measurement of charm flow with the STAR Heavy Flavor Tracker.", 
J. Phys. G: Nucl. Pari. Phys., 32, S571-S574, 2006. 

• Y. F. Zhang, S.Esumi, H.Huang, Y.Miake, S.Sakai and N.Xu, "Identifying bot- 
tom contribution in non-photonic electron spectra and V2 from Au-|-Au collisions at 
RHIC", Nucl. Phys. A, 783, 489-492, 2007. 

• H. D. Liu, Y. F. Zhang, C. Zhong and Z. Xu, "Extracting the charm cross-section 
from its semileptonic decay at RHIC." , Phys. Lett. B, 639, 441, 2006. 

• Y. F. Zhang, S. L. Huang, Z. R Zhang, J. Wu, "Study on "soft" and "hard" in- 
teractions in pp{pp) collisions using HIJING and PYTHIA." , Chinese Physics, 16, 
58-61, 2007. 

• Y. F. Zhang, S. L. Huang, Z. P. Zhang, J. Wu, "Simulation on "soft" and "hard" 
process in ppipp) collisions using HIJING and PYTHIA." , Journal of University of 
Science and Technology of China, 35(6), 821-824, 2005 (in Chinese). 

• Q. Shan, J. Wu, X. L. Wang, H. F. Chen, Y. F. Zhang, Z. P. Zhang, M. Shao, 
"Upper limit of the yield of di-omega in central Au-|-Au coUision at ^/{sNN) — 200 
GeV with HIJING.", High Energy Phys. and Nucl. Phys., 29, 1146-1149, 2005. 

• J. Adams et al, STAR Collaboration. "Pion, Kaon, Proton and Anti- proton Trans- 
verse Momentum Distributions from p -|- p and d -|- Au Collisions at y/s^^= 200 
GeV.", Phys. Lett. B, 616 8, 2005. 

• J. Adams et al., STAR Collaboration. "Hadronization Geometry and Charge- 
dependent Number Autocorrelations on Axial Momentum Space in Au-|-Au Col- 
hsions at ^/s^^ 130 GeV. " , Phys. Lett. B, 634 347, 2006. 



103 



J. Adams et al., STAR Collaboration. "Azimuthal Anisotropy in Au+Au Collisions 
at 200 GeV. ", Phys. Rev. C, 72 014904, 2005. 

J. Adams et ai, STAR Collaboration. "Minijet Deformation and Charge-independent 
Angular Correlations on Momentum Subspace (77, (p) in Au+Au Collisions at -^/s^— 
130 GeV. ", Phys. Rev. C, 73 064907, 2006. 

J. Adams et al, STAR Collaboration. "Experimental and Theoretical Challenges 
in the Search for the Quark Gluon Plasma: The STAR Collaboration's Critical 
Assessment of the Evidence from RHIC Collisions. ", Nucl. Phys. A, 757 102, 
2005. 

J. Adams et al, STAR Collaboration. "Multiplicity and Pseudorapidity Distribu- 
tions of Photons in Au-|-Au Collisions at y^s^= 62.4 GeV. ", Phys. Rev. Lett, 95 
062301, 2005. 

J. Adams et al, STAR Collaboration. "Multi-strange Baryon Elliptic Flow in 
Au+Au collisions at 200 GeV. ", Phys. Rev. Lett, 95 122301, 2005. 

J. Adams et al., STAR Collaboration. "Incident Energy Dependence of pt Correla- 
tions at RHIC. ", Phys. Rev. C, 72 044902, 2005. 

J. Adams et al., STAR Collaboration. "Transverse- momentum Correlations on 
(77, 0) from mean pt Fluctuations in Au+Au Collisions at -yjs^— 200 GeV. ", J. 
Phys. G, 32 L37, 2006. 

J. Adams et al., STAR Collaboration. "Directed Flow in Au+Au Collisions at 
^/s^^ 62 GeV. ", Phys. Rev. C, 73 034903, 2006. 

J. Adams et al., STAR Collaboration. "Proton- Lambda Correlations in Central 
Au+Au Collisions at ^/s^^ 200 GeV. " , Phys. Rev. C, 74 064906, 2006. 

J. Adams et al., STAR Collaboration. "Multiplicity and Pseudorapidity Distribu- 
tions of Charged Particles and Photons at Forward Pseudorapidity in Au+Au Col- 
hsions at y/s^^^ 62.4 GeV. ", Phys. Rev. C, 73 034906, 2006. 



104 



J. Adams et al., STAR Collaboration. "Strangelet Search at RHIC. ", arXiv: nucl- 
ex/0511047. 

M. Calderon de la Barca Sanchez et al., STAR Collaboration. "Open Charm Pro- 
duction from d + Au CoUisions in STAR. " , Eur. Phys. J. C, 43 187, 2005. 

A. A. P. Suaidc et al., STAR Collaboration. "Charm Production in the STAR 
Experiment at RHIC. " , Eur. Phys. J. C, 43 193, 2005. 

C. A. Gagliardi et al., STAR Collaboration. "Recent high-p-r Results from STAR. 
", Eur. Phys. J. C, 43 263, 2005. 

J. Adams et al., STAR Collaboration. "Identified Hadron Spectra at Large Trans- 
verse Momentum inp -\- p and d + Au CoUisions at ■y/s^= 200 GeV. " , Phys. Lett. 

B, 637 161, 2006. 

J. Adams et al., STAR Collaboration. "Measurements of Identified Particles at In- 
termediate Transverse Momentum in the STAR Experiment from Au-|-Au Collisions 
at ,/s^^2m GeV. ", arXiv: nucl-ex/0601042 . 

J. Adams et al., STAR Collaboration. "Forward Neutral Pion Production in p + p 
and d + Au Collisions at y/s^^^ 200 GeV. ", Phys. Rev. Lett, 97 152302, 2006. 

J. Adams et al., STAR Collaboration. "Direct Observation of Dijets in Central 
Au+Au Collisions at ^/s^^^ 200 GeV. ", Phys. Rev. Lett, 97 162301, 2006. 

J. Adams et al., STAR Collaboration. "Strange Baryon Resonance Production in 
200 GeVp + p and Au+Au CoUisions. ", Phys. Rev. Lett, 97 132301, 2006. 



J. Adams et al., STAR Collaboration. "The Energy Dependence of pt Angular 
Correlations Inferred from mean pt Fluctuation Scale Dependence in Heavy Ion 
CoUisions at the SPS and RHIC. ", J. Phys. G, 33 451, 2007. 

J. Adams et al., STAR Collaboration. "Identified Baryon and Meson Distributions 
at Large Transverse Momenta from Au-|-Au CoUisions at a/s^— 200 GeV. " , Phys. 
Rev. Lett, 97 152301, 2006. 



105 



J. Adams et ai, STAR Collaboration. "Scaling Properties of Hyperon Production 
in Au+Au Collisions at y/s^^ 200 GeV. ", arXiv: nucl-ex/0606014. 

J. Adams et a/., STAR Collaboration. "The Multiplicity Dependence of Inclusive 
Pt Spectra from p + p Collisions at 200 GeV. ", Phys. Rev. D, 74 032006, 



J. Adams et al, STAR Collaboration. "Delta phi Delta eta Correlations in Central 
Au+Au Collisions at 200 GeV. ", Phys. Rev. C, 75 034901, 2007. 



B. I. Abelev et al, STAR Collaboration. "Transverse Momentum and Centrahty 
Dependence of high-py Non-photonic Electron Suppression in Au+Au Collisions at 
200 GeV. ", arXtv: nucl-ex/0607012. 

B. I. Abelev et al., STAR Collaboration. "Strange Particle Production in p + p 
CoUisions at ./s^= 200 GeV. ", arXiv: nud-ex/0607033. 



B. I. Abelev et al., STAR Collaboration. "Neutral Kaon Interferometry in Au+Au Col- 
lisions at 200 GeV. ", Phys. Rev. C, 74 054902, 2006. 



B. 1. Abelev et al, STAR Collaboration. "Longitudinal Double-spin Asymmetry 
and Cross Section for Inclusive Jet Production in Polarized Proton Collisions at 
200 GeV. ", Phys. Rev. Lett, 97 252001, 2006. 

B. I. Abelev et al., STAR Collaboration. "Rapidity and Species Dependence of 
Particle Production at Large Transverse Momentum for d + Au CoUisions at ■Js^— 



200 GeV. " , arXiv: nucl-ex/0609021. 

B. I. Abelev et al., STAR Collaboration. "Partonic Flow and Phi-meson Production 
in Au + Au CoUisions at 200 GeV. " , arXiv: nucl-ex/0703033. 

B. I. Abelev et al., STAR Collaboration. "Charged Particle Distributions and Nu- 
clear Modification at High Rapidities in d + Au Collisions at ^Js^— 200 GeV. ", 



2006. 



arXiv: nucl-ex/0703033. 



106 



• B. I. Abelev et al., STAR Collaboration. "Mass, Quark- number, and -s/s^ De- 
pendence of the Second and Fourth Flow Harmonics in Ultra-Relativistic Nucleus- 
Nucleus Collisions. ", arXiv: nucl-ex/0701010. 



107 



References 



[AbrOl] 
[ABR03] 

[Acc04] 

[Aco03] 

[Ada03a] 

[AdaOSb] 

[Ada03c] 

[AdaOSd] 

[Ada04a] 

[Ada04b] 

[Ada04c] 

[Ada04d] 

[Ada04e] 

[Ada04f] 

[Ada04g] 

[Ada05] 

[Adc02] 

[Adc04] 

[Adl02a] 

[Adl02b] 

[AdlOSa] 

[Adl03b] 

[Adl04a] 

[Adl04b] 



M.C. Abreu et al. Phys. Lett. B, 521:195, 2001. 

A. Andronic, P. Braun-Munzinger, K. Redlich, and J. Stachel. Phys. Lett. B, 
571:36, 2003. 

A. Accardi and M. Gyulassy. J. Phys. G, 30:S969, 2004; Phys. Lett. B, 
586:244, 2004. 

D. Acosta et al. Phys. Rev. Lett., 91:241804, 2003. 
D. Adamova et al. Phys. Rev. Lett, 91:042301, 2003. 



J. Adams et a' 
J. Adams et a. 
J. Adams et a 
J. Adams et a' 
J. Adams et a. 
J. Adams et a 
J. Adams et a. 
J. Adams et a 
J. Adams et a 
J. Adams et a' 
J. Adams et a 
K. Adcox et a] 
K. Adcox et a 



Phys. Rev. Lett, 91:072304, 2003. 
arXiv: nucl-ex/ 03090 12, 2003. 
Phys. Rev. Lett, 91:172302, 2003. 
Phys. Rev. Lett, 92:112301, 2004. 
Phys. Rev. Lett, 92:052302, 2004. 
arXiv: nucl-ex/0406003, 2004. 
Phys. Rev. C, 70:044902, 2004. 
Phys. Rev. Lett, 93:252301, 2004. 
Phys. Rev. Lett, 93:252301, 2004. 
Phys. Rev. Lett, 92:171801, 2004. 
"STAR white paper." arXiv: nucl-ex/0501009, 2005. 
Phys. Rev. Lett, 88:192303, 2002. 

"PHENIX white paper." arXiv: nucl-ex/04 1 0003, 2004. 
C. Adler et al. Phys. Rev. Lett, 89:202301, 2002. 
C. Adler et al. Phys. Rev. C, 66, 2002. 
S.S. Adler et al. Phys. Rev. Lett, 91:182301, 2003. 
S.S. Adler et al. Phys. Rev. Lett, 91:241803, 2003. 
S.S. Adler et al. arXiv: nucl-ex/04 09015, 2004. 
S.S. Adler et al. Phys. Rev. Lett, 92:051802, 2004. 



108 



[Adl04c] S.S. Adler et al. Phys. Rev. C, 69:014901, 2004. 

[Adl05] S.S. Adler et al. Phys. Rev. Lett., 94:082301, 2005. 

[Afa02] S.V. Afanasiev et al. Phys. Rev. C, 66:054902, 2002. 

[AlaOl] J. Alam et al. Phys. Rev. C, 63:021901, 2001. 

[Alp02] S. Bethke. Nucl. Phys. Proc. Suppl, 121:74-81, 2003. 

[Altos] C. Alt et al. Phys. Rev. C, 68:034903, 2003. 

[Alv96] G.A. Alves et al. Phys. Rev. Lett., 77:2388, 1996. 

[Amm88] R. Ammar et al. Phys. Rev. Lett., 61:2185, 1988. 

[And03] M. Anderson et al. Nucl. Instr. Meth. A, 499:659, 2003. 

[Arm05] N. Armesto et al. Phys. Rev. D, 71:054027, 2005. 

[Arm06] N. Armesto et al. Phys. Lett. B, 637:362, 2006. 

[B76] F.W. Biisser et al. Nucl. Phys. B, 113:189, 1976. 

[Bar88] S. Barlag et al. Z. Phys. C, 39:451, 1988. 

[BCH03] E.L. Bratkovskaya, W. Gassing, and Stocker H. Phys. Rev. C, 67:054905, 2003. 

[BCS04] E.L. Bratkovskaya, W. Gassing, H. Stocker, and N. Xu. arXiv: nucl- 
th/0409047, 2004. 

[Bea04] I.G. Bearden et al. Phys. Rev. Lett., 93:102301, 2004. 

[Bic06] H. Bichsel. Nucl. Instr. Meth. A, 562:154, 2006. 

[Bjo83] J.D. Bjorken. Phys. Rev. D, 27:140, 1983. 

[BKG03] S. Batsouli, S. Kelly, M. Gyulassy, and J.L. Nagle. Phys. Lett. B, 557:26, 
2003. 

[Bot90] O. Botner et al. Phys. Lett. B, 1990. 

[BPLF04] G. G. Barnafoldi, G. Papp, P. Levai and G. Fai. J. Phys. G, 30:81125, 2004. 

[BraOl] P. Braun-Munzinger. Nucl. Phys. A, 681:119c, 2001. 

[Bus04] P.J. Bussey et al. arXiv: hep-ex/0408020, 2004. 

[Gat04] E. Gattaruzza and D. Treleani. Phys. Rev. D, 69:094006, 2004. 

[Gss89] J.G. GoUins, D.E. Soper and G. Sterman, in Perturbative QCD, edited by 
A.H. Mueller (World Scientific, Singapore, 1989), p.l. 



109 



[Col98] J.C. Collins. Phys. Rev. D, 66:094002, 1998. 
[Cas04] J. CastilUo et al. J. Phys. G, 30:81207, 2004. 

[CER] CERN. "littp://wwwasd.wcb.cern.ch/wwwasd/cernlib/.". 

[CLEO04] CLEO Collaboration, J. Yclton et al. Presented at Heavy Quarks and Leptons, 
San Juan, Puerto Rico, June 4, 2004 

[CNV05] M. Cacciari, P. Nason and R. Vogt. Phys. Rev. Lett., 95:122001, 2005. 

[CQM05] H. B. Zhang et al. Quark Matter 2005, Budapest, Hungary, 4-9 Aug. 2005, 
arXiv: nucl-ex/05 10063; J. Bielcik et al. arXiv: nucl-ex/0511005. 

[Cro73] J.W. Cronin et al. Phys. Rev. Lett, 31:1426, 1973. 

[DDZ06] X. Zhu et al. arXiv: hep-ph/0604178. 

[DES04] X. Dong, S. Esumi, P. Sorensen, N. Xu, and Z. Xu. Phys. Lett. B, 597:328, 
2004. 

[DFS03] D. de Florian and R. Sassot. Phys. Rev. D, 69:074028, 2004. 

[DG03] M. Djordjevic and M. Gyulassy. Phys. Lett. B, 560:37, 2003. 

[DG04a] M. Djordevic and M. Gyulassy. J. Phys. G, 30:81183, 2004. 

[DC04b] M. Djordjevic and M. Gyulassy. Nucl. Phys. A, 733:265, 2004. 

[DGV06] M. Djordjevic, M. Gyulassy, R. Vogt and S. Wicks. Phys. Lett. B, 632:81, 
2006. 

[DGW05] M. Djordjevic, M. Gyulassy, and S. Wicks. Phys. Rev. Lett, 94:112301, 2005. 

[DKOl] Y.L. Dokshitzer and D.E. Kharzeev. Phys. Lett. B, 519:199, 2001. 

[DKS03] G. Dissertori, 1. Knowlcs, and M. Schmelling. Quantum Chromodynamics - 
High Energy Experiments and Theory. Oxford Universidy Press, 2003. 

[EKS98] K. J. Eskola, V. J. Kolhinen and C. A. Salgado. Eur. Phys. J. C, 9:61, 1999. 

[FK02] Z. Fodor and S.D. Katz. J. High Energy Phys., 0203:014, 2002. 

[GlauOl] P. Shukla arXiv: nucl-th/01 12039, 2001. 

[FMN97] S. Frixione, M.L. Mangano, P. Nason, and G. Ridolfi. arXiv: hep-ph/9702287, 
1997. 

[HouOl] P. Houvinen et al. Phys. Lett B, 503:58, 2001. 

[GK04] V. Greco and CM. Ko. Phys. Rev. C, 70:024901, 2004. 



110 



[GKL02] M.I. Gorenstcin, A. P. Kostyuk, McLerran L., H. Stocker, and W. Greiner. J. 
Phys. G, 28:2151, 2002. 

[GKR04] V. Greco, C.M. Ko, and R. Rapp. Phys. Lett B, 595:202, 2004. 

[GROl] L. Grandchamp and R. Rapp. Phys. Lett. B, 523:60, 2001. 

[Gra04] L. Grandchamp. Private communications, 2004. 

[Gup98] R. Gupta. arXiv: hep-lat/9807028, 1998. 

[Hbz03] H. B. Zhang. Ph.D. Thesis, Yale University, 2003. 

[HY04] R.C. Hwa and C.B Yang. Phys. Rev. Lett, 93:082302, 2004. 

[HYP06] J. Adams et al. arXiv: nucl-ex/0606014. 

[HGR06] H. Hees, V. Greco and R. Rapp. Phys. Rev. C, 73:034913, 2006. 

[Joh02] I. Johnson. Ph.D. Thesis, U.C. Davis, 2002. 

[Kar02a] F. Karsch. Lect Notes Phys., 583:209, 2002. 

[Kar02b] F. Karsch. Nucl. Phys. A, 698:199, 2002. 

[KH03] P.F. Kolb and U. Heinz. arXiv: nucl-th/0305084, 2003. 

[KLOl] D. Kharzeev and E Levin. Phys. Lett B, 523:79, 2001. 

[KLPOO] F. Karsch, E. Laermann, and A. Peikert. Phys. Lett B, 478:447, 2000. 

[Kod91] K. Kodama et al. Phys. Lett B, 263:573, 1991. 

[Kop03] B. Z. Kopeliovich. Phys. Rev. C, 68:044906, 2003. 

[Kos02] Kostyuk et al. Phys. Lett B, 524:265, 2002. 

[LG96] Z. Lin and M. Gyulassy. Phys. Rev. Lett, 77:1222, 1996. 

[Lin95] Z. Lin et al. Phys. Rev. C, 51:2177, 1995. 

[Lin96] Z. Lin. Ph.D. Thesis, Columbia University, 1996. 

[LXY04] X. Y. Lin. arXiv: hep-ph/0412124. 

[LXY06] X. Y. Lin. arXiv: hep-ph/ 0602067. 

[LK02] Z. Lin and C.M. Ko. Phys. Rev. Lett, 89:202302, 2002. 

[Llo04] W.J. Llope et al. Nucl. Instr. Meth. A, 522:252, 2004. 

[LM03] Z. Lin and D. Molnar. Phys. Rev. C, 68:044901, 2003. 



Ill 



[Lqe06] J. Adams et al. Phys. Rev. Lett, 97:152301, 2006. 

[LZZX06] H. Liu, Y. Zhang, C. Zhong and Z. Xu. Phys. Lett. B, 639:441, 2006. 

[MNR93] M.L. Mangano, P. Nason, and G. Ridolfi. Nucl. Phys. B, 405:507, 1993. 

[MS86] T. Matsui and H. Satz. Phys. Lett. B, 178:416, 1986. 

[MT05] G.D. Moore and D. Teaney. Phys. Rev. C, 71:064904, 2005. 

[MV03] D. Molnar and S.A. Voloshin. Phys. Rev. Lett, 91:092301, 2003. 

[NLC03] G.C. Nayak, M.X. Liu, and F. Gooper. Phys. Rev. D, 68:034003, 2003. 

[PbmOO] P. Braun-Munzinger and J. Stachel. Phys. Lett B, 490:196, 2000. 

[PDG02] K. Hagiwara et al. Phys. Rev. D, 66:010001, 2002. 

[PDG04] S. Eidelman et al. "Review of Particle Physics." Phys. Lett B, 592:1+, 2004. 

[PDG06] W. M. Yao et al. J. Phys. G, 33:1, 2006. 

[PHEGAA] S. S. Adler et al. Phys. Rev. Lett, 94:082301, 2005. K. Adcox et al. Phys. 
Rev. Lett, 88:192303, 2002. 

[PHEGPP] S. S. Adler et al. Phys. Rev. Lett, 97:252002, 2006. 

[PHEGRAA] S. S. Adler et al. Phys. Rev. Lett, 96:032301, 2006. 

[PHEMu] S. S. Adler et al. arXiv: hep-ex/0609032. 

[PhpiO] S. S. Adler et al. Phys. Rev. Lett, 91:072301, 2003. 

[Pic95] A. Pich. arXiv: hep-ph/9505831, 1995. 

[Polem] G. S. Huang ct al. Phys. Rev. Lett, 94:011802, 2005; E. M. Aitala et al. Phys. 
Rev. Lett, 80:1393, 1998. 

[PV98] A.M. Poskanzer and S.A. Voloshin. Phys. Rev. C, 58:1671, 1998. 

[RapOl] R. Rapp. Phys. Rev. C, 63:054907, 2001. 

[Rap06] R. Rapp and H. Hccs. J. Phys. G, 32:S351, 2006. 

[RP03] J. Raufeisen and J.-G. Peng. Phys. Rev. D, 67:054008, 2003. 

[RRL03] I.V. Rakobolskaya, T.M. Roganova, and Sveshnikova L.G. Nucl. Phys. B (Proc. 
Suppl), 122:353, 2003. 

[RS03] R. Rapp and E.V. Shuryak. Phys. Rev. D, 67:074036, 2003. 

[Rua04a] L.J. Ruan et al. J. Phys. G, 30:81197, 2004. 



112 



[Rua04b] L.J. Ruan. Ph.D. Thesis, University of Science and Technology of China, 2004. 

[Se04] B. MuUer arXiv: nucl-th/0404015. 

[Sha02] M. Shao et al. Nucl. Instr. Meth. A, 492:344, 2002. 

[Sha05] M. Shao et al. J. Phys. G, 31:S85, 2005. 

[Sha06] M. Shao et al. Nucl. Instr. Meth. A, 558:419, 2006. 

[SjOl] T. Sjostrand et al. Computer Physics Commmun., 135:238, 2001. 

[Sor03] P.R. Sorensen. Ph.D. Thesis, UC Los Angeles, 2003. arXiv: nucl-ex/0309003. 

[Sor04a] P.R. Sorensen. Private communications., 2004. 

[Sor04b] P.R. Sorensen. STAR internal talk, 2004. 

[Sor99] H. Sorge. Phys. Rev. Lett, 82:2048, 1999. 

[SSH93] E. Schnedermann, J. SoUfrank, and U. Heinz. Phys. Rev. C, 48:2462, 1993. 

[SS07] S. Sakai. Quark Matter 2006 proceeding, Shanghai, China. 

[STA03] K.H. Ackermann et al. Nucl. Instr. Meth. A, 499:624, 2003. 

[STA04] STAR TOF Collaboration. STAR TOF Proposal, 2004. 

[STA05] STAR HFT Collaboration. STAR HFT Proposal, 2005; Z. Xu et al. J. Phys. G: 
Nucl. Part. Phys. :32, S571, 2006; A. Rose et al. Quark Matter 2006 proceeding, 
Shanghai, China. 

[STAc05] J. Adams et al. Phys. Rev. Lett., 94:062301, 2005. 

[STACRAA] J. Adams et al. arXiv: nucl-ex/0607012. 

[Sto04] R. Stock. J. Phys. G, 30:S633, 2004. 

[Tai04] A. Tai et al. J. Phys. G, 30:S809, 2004. 

[Tav87] S.P.K. Tavernier. Rep. Prog. Phys., 50:1439-1489, 1987. 

[Tof05] J. Adams et al. Phys. Lett. B, 616:8, 2005. 

[TSROl] R.L. Thews, M. Schroedter, and J. Rafelski. Phys. Rev. C, 63:054905, 2001. 

[Vit03] I Vitev. Phys. Lett. B, 562:36, 2003. 

[Vit06] A. Adil and 1. Vitev. arXiv: hep-ph/ 061 1109. 

[Vog02] R. Vogt. arXiv: hep-ph/0203151, 2002. 

[Vog03] R. Vogt. Int. J. Mod. Phys. E, 12:211, 2003. 



113 



[Vog04] R. Vogt. Private communications, 2004. 

[Vogt04] R. Vogt. Phys. Rev. C, 70:064902, 2004. 

[WanOO] X.N. Wang. Phys. Rev. C, 61:064910, 2000. 

[Wan04] F. Wang et al. J. Phys. G, 30:81299, 2004. 

[WG91] X.N. Wang and M. Gyulassy. Phys. Rev. D, 44:3501, 1991. 

[Wic05] S. Wicks et al. arXiv: nucl-th/05 12076. 

[Won94] C.-Y. Wong. Introduction to High-Energy Heavy-Ion Collisions. World Scien- 
tific Publishing Co. Pte. Ltd., 1994. 

[Xin05] X. Dong. Ph.D. Thesis, USTC Hefei, China, 2005. arXiv: nucl-ex/0509011. 

[Xu03] N. Xu and Z. Xu. Nucl. Phys. A, 715:587c, 2003. 

[Xu04] Z. Xu et al. arXiv: nucl-ex/OJ^llOOl, 2004. 

[Yan07] Y. Lu. Ph.D. Thesis, CCNU Wuhan, China, 2007. 

[YF06] Y. F. Zhang. J. Phys. G: Nucl. Part. Phys., 32:8529-8532, 2006. 

[YFHP06] Y. F. Zhang, 8. Esumi, H. Huang, Y. Miake, 8. 8akai and N. Xu. Nucl. Phys. 
A, 783:489, 2007 

[ZCK05] B. Zhang, L.W. Chen, and CM. Ko. Phys. Rev. C, 72:024906, 2005. 
[ZGK99] B. Zhang, M. Gyulassy, and CM. Ko. Phys. Lett. B, 455:45, 1999. 



114 



